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CHAPTER  I 


INTRODUCTION 


1 . 1  BACKGROUND 

The  rapid  growth  in  air  transportation  over  the  last  four 
decades  and  the  increases  in  payload  of  commercial  aircraft  have 
led  to  severe  deterioration  of  airport  pavements.  Invariably, 
these  pavements  are  rehabilitated  to  preserve  or  increase  load 
carrying  capabilities  or  to  maintain  a  smooth  riding  surface. 
Common  practice  is  to  overlay  the  cracked  pavement  with  a  matting 
of  asphaltic  material,  or  in  some  instances,  with  portland  cement 
concrete.  In  any  case,  a  new  problem  usually  develops:  the  new 
overlay,  after  a  short  time,  begins  to  crack  in  the  same  pattern 
as  the  pavement  that  was  overlaid;  this  constitutes  the 
phenomenon  of  reflection  cracking. 

Large  sums  of  money  are  spent  each  year  by  pavement 
management  groups  to  seal  these  cracks  against  the  ingress  of 
water  and  the  growth  of  vegetation;  other  pavements  that  have 
been  allowed  to  go  unsealed  have  in  a  very  short  time 
necessitated  complete  reconstruction  efforts.  While  many 
remedial  measures  have  been  taken  to  prevent  or  retard  reflection 
cracking,  no  generally  applicable  method  has  been  found  that 
would  apply  to  all  pavement  types.  One  of  the  most  promising  of 
those  measures  has  emerged  to  be  the  addition  of  an  engineering  type 
of  fabric  to  the  cracked  pavement  before  application  of  the 
bituminous  overlay.  Although  many  millions  of  dollars  are 
expended  annually  in  the  certain  belief  that  installation  of 
fabric  will  retard  cracking  of  airport  pavements,  much  of  this 
sum  is  wasted  because  of  failure  of  the  fabric  to  perform  due  to 
the  lack  of  adequate  construction  controls  and  design  criteria. 


1.2  OBJECTIVE  AND  SCOPE 

The  objective  of  this  research  effort  is  the  development  of 
a  mechanistic  design  methodology  that  permits  the  effective  use 
of  engineering  fabrics  on  pavement  overlay  systems  for  the 
prevention  of  reflection  cracking. 

The  scope  of  this  research  project  includes  literature 
search  and  review  of  analytical  models  related  to  the  use  of 
fabrics  as  crack  arresting  systems.  Based  on  this  review, 
functional  criteria  are  established  and  existing  analytical 
procedures  modified  or  improved  for  the  analysis  and  design  of 
these  systems.  Prior  results  on  the  use  of  fabrics  on  pavements 
are  augmented  by  additional  laboratory  work  and  theoretical 
studies.  A  modular  computer  program  is  also  developed  and 
acceptance  tested  and  documented  in  the  form  of  a  User's  Manual 


and  an  Operations  Manual.  The  scope  of  this  development  effort 
also  includes  sufficient  laboratory  testing  of  small  scale  models 
for  representative  types  of  pavements  with  varying 
characteristics  of  engineering  fabrics  included.  Field 

verification  of  the  structural  model  will  be  performed  before  _a 
design  manual  is  formulated  from  the  total  results  of  this 
effort. 

1.3  OUTLINE  OF  STUDY 

This  study  is  divided  into  two  volumes.  Volume  I  describes 
the  development  of  the  design  method  including  an  appendix 
containing  the  user  guide  which  also  includes  guidelines  for 
development  of  the  input  parameter  required  by  the  method. 

Volume  I  is  divided  into  the  following  sections: 

CHAPTER  II:  Review  of  Literature 

CHAPTER  III:  Mechanistic  Models  for  Geotextile  Systems 
CHAPTER  IV:  Design  Methods 
CHAPTER  V:  Laboratory  Testing 

CHAPTER  VI:  Data  Analysis 

CHAPTER  VII:  Conclusion  and  Recommendations 
APPENDIX:  User  Guide 

Volume  I  documents  the  research  efforts  of  the  first  phase 
of  this  study. 

Volume  II  will  include  activities  of  Phase  III,  "Field 
Verification",  which  is  currently  underway.  The  design  model 
presented  in  this  report  may  require  modification  and/or 
calibration  as  a  result  of  the  field  studies.  It  should  be  noted 
that  since  the  field  phase  is  still  under  investigation.  Volume  I 
is  limited  at  this  stage  to  a  computer  program  input  guide  rather 
than  a  full  design  model.  The  final  design  model  will  be 
presented  at  the  conclusion  of  all  activities. 


CHAPTER  II 


REVIEW  OF  LITERATURE 


2.1.  INTRODUCTION 

Reflection  cracking  is  the  cracking  of  a  resurface  or 
overlay  above  underlying  cracks  or  joints.  This  cracking  occurs 
in  overlays  of  both  flexible  and  rigid  pavements  and  is  a  major 
cause  of  future  pavement  distress  including  spalling,  surface 
water  infiltration  to  underlying  base  and  subgrade  layers,  and  a 
general  reduction  in  the  stiffness  of  the  pavement  structure. 
Reflective  cracks  require  labor  intensive  operations  for  crack 
sealing  and  patching,  thus  becoming  a  significant  maintenance 
expense  item. 

The  problem  of  reflection  cracking  is  not  a  new  one  to  the 
pavement  engineer.  Since  the  early  1950's  many  different 
materials ,  methods  and  techniques  have  been  tried  to  prevent  or 
at  least  delay  reflection  cracking.  Most  of  these  efforts  have 
been  concerned  with  asphalt  concrete  overlays  over  existing 
Portland  Cement  Concrete  pavements,  where  existing  cracks  or 
joints  are  usually  reflected  through  the  asphalt  overlay  within  a 
1  year  period  (1).  Early  research  recognized  that  the  probable 
cause  of  reflection  cracking  was  movement  of  some  form  in  the 
underlying  pavement  at  existing  cracks  and  joints.  This  movement 
can  result  from  both  traffic  and  environmentally  induced  forces. 
The  movement  includes  differential  vertical  movement  and  thermal 
or  moisture  induced  expansion,  contraction  or  distortion 
(curling)  at  underlying  joints  and  cracks.  Because  the  overlay 
is  bonded  to  the  existing  pavement,  movement  at  underlying  joints 
or  cracks  induces  stresses  in  the  overlay.  If  sufficiently  high, 
these  stresses  cause  fracturing  or  cracking  of  the  overlay.  If 
the  induced  stresses  do  not  exceed  the  yield  strength  of  the 
overlay  material,  cracking  could  still  develop  as  the  result  of 
cyclic  load  applications  which  produce  fatigue  fracturing  of  the 
asphalt  concrete.  Techniques  or  measures  which  have  been  used  to 
delay  reflection  cracking  are:  (i)  Bond  Breakers;  (ii) 
Cushions;  (iii)  Rubber  -  Asphalt  Interlayer  (SAMI);  (iv) 
Fabrics;  (v)  Modifying  Existing  Pavement;  (vi)  Stronger 
Overlays . 

Currently  the  technique  subjected  to  the  most  extensive 
field  testing  to  determine  its  effectiveness  in  delaying 
reflection  cracking  is  the  placement  of  engineering  fabrics  over 
cracks  or  joints  prior  to  overlay.  There  are  numerous 
experimental  projects  of  this  type  throughout  the  country. 
Reference  (1)  contains  a  tabulation  of  experimental  projects 
using  "Petromat"  fabric  including  location,  test  features, 
evaluation  procedures,  and  observations  to  date.  "Petromat"  is  a 
nonwoven  polypropylene  fabric  produced  by  Phillips  Petroleum  and 
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has  been  the  most  widely  used  fabric  to  date.  There  are  at  least 
6  other  major  U.S.  Corporations  manufacturing  fabrics  from  such 
materials  as  nylon,  polyester,  polypropylene,  polyvinylidene 
chloride,  or  fiberglass. 

Fabrics  are  available  in  rolls  and  are  placed  as  sheets  tack 
coated  to  the  existing  pavement  prior  to  overlay.  Recently 

manufacturers  have  started  producing  2  to  3  foot  (3/4  to  Ira)  wide 
rolls  so  that  the  fabric  can  be  placed  as  strips  over  joints  or 
cracks.  These  narrow  rolls  sometimes  have  an  adhesive  (usually 
rubberized  asphalt)  on  one  side  so  they  can  be  placed  without 
tack  coating  the  existing  surface.  The  strip  type  fabrics  are 
usually  used  on  rigid  pavements  while  sheet  fabrics  which  cover 
the  entire  area  are  used  on  flexible  pavements.  In  a  Virginia 
project  3  foot  (0.91m)  wide  strips  of  non-woven  polypropylene 
fabrics  were  placed  over  a  composite  pavement  at  existing 

transverse  reflection  cracks  (2).  Joint  spacing  of  the  PCC  base 
was  30  feet  (9.1m).  An  emulsion  tack  coat  was  applied  under  the 
fabric  and  overlaid  with  1  1/4  inches  (32mm)  of  asphalt  concrete. 
After  three  months  under  traffic  many  of  the  joints  were 
reflected  through  the  second  overlay  although  there  was  somewhat 
more  cracking  in  an  adjacent  section  where  no  fabric  has  been 
used.  Virginia  constructed  other  projects  with  fabric  placed 
directly  on  PCC  slabs  at  joints  and  cracks  prior  to  overlay  and 
attempted  to  correlate  the  fabric's  crack  prevention  performance 
with  differential  vertical  movement  or  load  transfer  under  18  kip 
axle  loads  at  the  joints.  The  Virginia  study  found  a  significant 
correlation  between  fabric  performance  and  joint  load  transfer. 
The  researchers  stated  that  "...  it  is  likely  that  many  of  the 
20  fabric-treated  joints  that  were  uncracked  and  had  a 

differential  deflection  of  0  (100%  load  transfer)  were 

working  joints  where  the  fabric  served  its  intended  purpose  of 
reducing  overlay  stresses  to  the  point  that  no  cracking  occurred. 
Conversely,  it  is  likely  that,  for  those  joints  that  had  higher 
differential  deflections,  the  fabric,  a  thin  sheet,  could  not 
sufficiently  distribute  the  shear  stresses  and  was  thus  unable  to 
reduce  reflection  cracking  significantly."  The  study  also 

concluded  that  at  differential  deflections  greater  than  .002  in 
(.05mm)  reflection  cracks  form  very  early.  Lower  differential 
deflections  delay  cracking,  but  cracks  will  occur  as  the 
magnitude  and  frequency  of  wheel  loads  increase. 

The  Virginia  study  appears  to  validate  the  fatigue  fracture 
mechanism  for  reflection  cracking  under  traffic  load  first 
presented  by  Majidzadeh,  et  al,  (3).  The  Virginia  results  are 
very  important  and  may  explain  why  other  states  have  had 
difficulty  in  assessing  fabric  reflection  crack  prevention 
performance.  A  recent  State  of  Art  review  conducted  by  FAA 
summarizing  highway  experience  with  fabrics  stated  that  no 
definite  conclusions  could  be  drawn  from  the  collective  results 
(4).  In  some  instances  apparently  the  fabric  (Petromat  and 
others)  appeared  to  be  effective,  but  in  apparently  similar 
situations  at  other  locations,  results  would  contradict  previous 
conclusions.  Based  upon  Army  Corps  of  Engineers  experience. 
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there  is  some  evidence  that  fabrics  are  more  effective  on  asphalt 
over  asphalt  than  on  asphalt  over  jointed  concrete  in  retarding 
reflection  cracking  (5;.  However,  fabric  placed  beneath  a  1  3/4 
in  (44mm)  overlay  over  an  existing  flexible  pavement  in  an 
Arizona  study  did  not  perform  favorably.  This  technique  was  not 
among  the  five  most  successful  methods  for  reducing  reflection 
cracking  at  the  Arizona  site  (6). 

The  reflection  cracking  studies  and  field  experimental 
projects  to  date  have  generally  been  of  an  empirical  nature  with 
little  control  or  even  identification  of  the  parameters  known  to 
affect  cracking.  Characterization  of  the  existing  pavement  in 
terms  of  joint  width,  load  transfer,  crack  spacing,  crack  and 
joint  opening  under  known  temperature  conditions,  and  deflection 
under  load  has  usually  not  been  part  of  these  studies. 
Obviously,  certain  crack  prevention  treatments  are  sensitive  to 
some  of  these  factors  as  shown  in  the  Virginia  study  where  fabric 
performance  was  related  to  load  transfer.  Unfortunately,  past 
research  has  not  established  the  quantitative  relationship 
between  these  factors  and  success  or  failure  of  the  preventative 
techniques . 


2.2  Evaluation  of  Existing  Reflection  Cracking  Models 

Within  the  last  10  years,  several  theoretical  (mathematical) 
models  have  been  developed  to  analyze  and  predict  the  occurrence 
of  reflection  cracking.  All  of  the  models  consider  the  same 
mechanisms  as  previously  noted  (e.g. ,  reflection  cracking  is 
caused  by  differential  horizontal  or  vertical  movements  in  the 
underlying  layer).  The  models  differ  in  the  methods  for 
predicting  the  magnitude  of  underlying  layer  movements,  on  the 
magnitude  of  stresses  induced  Into  the  overlay  by  the  movements, 
and  in  the  response  of  the  overlay  to  stress  state  (sudden 
fracture  vs.  fatigue  fracture).  The  following  discussions 
summarize  some  of  the  reflection  cracking  models  which  currently 
exist  and  present  an  evaluation  of  each  model's  limitations. 


2.2.1  Ohio  State  University  (OSU)  -  Ultimate  Strength  Model 

This  model,  developed  by  researchers  at  Ohio  State 
University  (OSU)  (7) ,  is  a  nomograph  procedure  for  predicting 
Asphalt  Conrete  (A.C.)  overlay  stresses  over  joints  or  cracks, 
resulting  from  thermally  induced  movements  in  underlying  portland 
cement  concrete  (PCC)  slabs.  Separate  stress  analyses  are 
performed  for  horizontal  slab  movements,  due  to  seasonal  changes 
in  average  slab  temperature,  and  vertical  slab  movements 
(curling)  which  occur  due  to  differential  vertical  temperature  of 
the  slabs.  Curling  is  the  response  state  where  the  top  of  the 
PCC  slab  is  colder  than  the  bottom  of  the  slab. 

The  horizontal  movement  of  the  PCC  slab  as  a  function  of  the 
change  in  slab  temperature  is  calculated  using  an  average  value 

for  the  friction  coefficient,  similar  to  the  calculation  for 
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determining  temperature  reinforcement  in  jointed  reinforced 
concrete  pavements.  This  model  neglects  the  resistance  to  joint 
movement  provided  by  the  uncracked  overlay  which  is  bonded  by  a 
tack  coat  to  the  underlying  slab.  The  OSU  model  assumes  that 
this  resistance  is  small  and  that  thin  overlays  do  not  affect  the 
movement  of  the  Joint  due  to  temperature  change.  The  joint 

dimension*  the  thickness,  and  the  modulus  of  elasticity  of  the 

overlay  and  the  slabs  are  input  into  a  finite  element  model  to 
determine  the  overlay  stresses. 

The  effect  of  vertical  movements  on  the  overlay  due  to  slab 
curling,  similar  to  the  horizontal  joint  movement,  is  also  based 
upon  the  premise  that  the  thin  overlays  do  not  affect  the  curling 
of  slabs  significantly.  Thus  the  restraint  against  curling  of 
the  slabs  provided  by  the  uncracked  overlay  is  again  neglected. 
This  important  assumption  permits  the  curved  shapes  of  slabs  to 
be  predicted  utilizing  a  computer  simulation  (program  called 

"PLATES")  of  the  Westergaard  solution  for  temperature 

differentials  between  top  and  bottom  of  the  slab.  Curling 

induced  overlay  stresses  are  estimated  on  the  assumption  that  the 
overlay  takes  the  slope  illustrated  in  Figure  1.  The  radius  of 
curvature  of  the  overlay,  R,  can  be  estimated  from  the  joint 
width, j, and  edge  slope,  e  ,  calculated  from  the  "PLATES"  program: 

J 

R  *  -  (1) 

28 

In  turn,  overlay  stresses  can  be  calculated  from  the  expression: 

Eov  hov 

<rov  =  (2) 

j 

where 

Eov  *  overlay  stiffness 
hov  =  overlay  thickness 

Equation  (2)  is  derived  from  the  basic  strength  of  materials  for 
pure  bending: 

«(u)  =  u/R 

«(u)  =  axial  strain  at  distance  u  from 
the  neutral  axis 

where 

R  =  radius  of  curvature 

Since  equation  (2)  is  derived  from  pure  bending,  symmetric 
bending  of  the  overlay  with  tension  at  the  top  and  compression  at 
the  bottom  is  implied. 


The  OSU-Ult <raate  Strength  Model  presents  an  easy  nomograph 
procedure  for  determining  overlay  stresses  from  thermal  movements 
of  underlying  PCC  slabs.  However,  the  accuracy  of  the  stress 
computation  is  suspect  for  the  following  reasons: 

a)  Restraint  imposed  by  the  uncracked  overlay  against  slab 
movements  (both  horizontal  and  curling)  is  not 
considered.  Thus,  the  calculated  force  in  the  overlay  at 
the  time  of  cracking  is  probably  incorrect. 

b)  The  overlay  stresses  due  to  horizontal  joint  movement 
appear  low  and  should  be  validated  by  additional  finite 
element  investigation. 

c)  The  tack  coat  bonding  stress  values  also  seem  low  and 
should  be  established  by  a  laboratory  investigation  which 
considers  temperature,  tack  coat  type  and  amount,  and 
roughness  of  the  PCC  slab. 

d)  The  Westergaard  analysis  used  to  predict  slab  curling 
neglects  the  weight ‘of  slab  and  overlay  which  would  tend 
to  reduce  the  curl. 

e)  The  simplified  analysis  of  overlay  stresses  due  to 

curling  should  be  verified  by  finite  element  analysis. 
The  fact  that  curling  introduces  a  horizontal  joint 
opening  is  neglected;  this  horizontal  movement  could  be 
significant  and  thus  change  the  stress  state 
considerably. 

f)  The  model  is  not  capable  of  assessing  the  effects  of 
crack  prevention  measures  upon  stresses  in  the  overlay. 

g)  The  model  does  not  present  recommendations  for  selecting 
design  parameters  such  as  seasonal  temperature  change  T, 
vertical  temperature  gradient  T,  asphalt  concrete 
modulus,  and  asphalt  concrete  strength. 

2.2.2  ARE  -  Ultimate  Strength  Model 

Austin  Research  Engineers  (8)  have  developed  a  procedure  for 
reflection  crack  stress  or  strain  analysis.  Two  different 
failure  modes  are  considered.  The  first  is  an  opening  mode 
(Figure  2)  due  to  horizontal  movements  of  the  underlying  PCC 
slab,  resulting  from  a  seasonal  temperature  change.  Both  joints 
or  cracks  without  steel  reinforcement,  or  cracks  with  steel 
reinforcement  (such  as  CRCP) ,  can  be  analyzed  for  horizontal 
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movement.  The  second  is  a  shearing  mode  (Figure  2)  resulting 
from  a  differential  deflection  across  the  joint  or  crack  as  the 
traffic  load  moves  across  the  discontinuity. 

A  number  of  assumptions  have  been  made  in  developing  the 
model,  including:  the  materials  are  elastic  in  response; 
temperature  variations  are  uniformly  distributed  in  the  existing 
concrete  slab  (no  curling);  concrete  movement  is  continuous  with 
slab  length;  and  movement  is  uniform  with  depth  in  a  particular 
layer. 

The  ARE-Ultiraate  Strength  Model  has  been  computerized 
(program  called  RFLCRl)  which  minimizes  difficulties  in  using  the 
model.  The  model  is  the  most  versatile  procedure  currently 
available  in  that  it  can  consider  slab  or  overlay  reinforcement, 
bond  breakers,  and  granular  cushions  (shear  failure  analysis 
only).  However,  the  simplifications  in  the  model,  which  permit 
strain  calculation  without  the  use  of  analytical  computations  of 
stress  distribution,  have  not  been  validated. 

Although  the  force  magnitudes  may  be  reasonable,  the  assumed 
simplified  distribution  of  stresses  within  the  overlay  for  both 
opening  and  shear  failure  modes  is  very  suspect  since  no 
concentration  of  stresses  at  the  joint  tip  is  considered.  Other 
less  significant  questions  regarding  the  ARE  model  include: 

a)  Characterization  of  the  existing  pavement  by  joint 
opening  measurements  over  a  certain  temperature  range 
cannot  necessarily  be  extrapolated  to  a  different  design 
temperature  range.  For  example,  restraint  exhibited 
between  70  and  50  degrees  F  (21  and  10  degrees  C)  may  not 
identify  the  restraint  between  70  and  20  degrees  F  (21 
and  -7  degrees  C) . 

b)  The  assumed  value  for  bonding  stress  between  overlay  and 
slab  is  very  important  to  the  anslysis  since  it 
establishes  the  gage  length  over  which  the  overlay  force 
at  the  joint  is  distributed.  The  suggested  values  need 
to  be  validated  experimentally. 

c)  The  concept  that  a  bond  breaker  reduces  overlay  strain  by 
merely  increasing  the  gage  length  for  force  transfer 
should  be  validated  by  analytical  investigation  of  stress 
distribution. 


d)  Load  transfer  is  determined  from  preoverlay  measurements 
without  any  adjustment  for  the  effect  of  the  overlay. 
This  effect  may  not  be  negligible.  Also  load  transfer  is 
probably  load  and  temperature  dependent. 


(e)  The  temperature  for  determination  of  dynamic  modulus  in 
the  shear  model  is  not  specified.  By  the  model,  a  high 
temperature  would  be  critical,  since  larger  strains  would 
result.  However,  the  allowable  strain  is  likely  to  be 
temperature  dependent. 


2.2.3  Ohio  State  University  (OSU)  -  Fracture  Mechanics  Model 

Fracture  mechanics  has  been  utilized  to  develop  a  reflection 
crack  propagation  model  for  asphalt  overlays  over  PCC  slabs  (3, 
9) .  The  model  considers  only  traffic  induced  fatigue  cracking 
resulting  from  differential  deflection  at  slab  joints  or  cracks. 

The  first  step  in  applying  fracture  mechanics  principles  was 
to  identify  the  fracture  mode  or  modes  associated  with  crack 
initiation  and  extension  (Figure  2).  A  finite  element  analysis 
of  full  scale  pavements  predictd  the  asphaltic  concrete  overlay 
to  be  in  compression,  thus  leading  to  the  conclusion  that  the 
opening  mode  (Mode  I)  type  fracture  does  not  occur. 
Additionally,  the  computer  analysis  predicted  that  there  would  be 
significant  relative  vertical  displacement  (Mode  II)  across  the 
joint  when  the  load  edge  is  placed  over  the  joint.  These 
conclusions  led  to  the  hypothesis  that  load  induced  reflection 
cracking  is  a  result  of  general  or  mixed  mode  fracture  of  the 
bituminous  material  occurring  under  the  simultaneous  interaction 
of  -K1  (negative  ^x) ,  K2  and  K3.  Laboratory  testing  of  2  and  3 
dimensional  model  overlay  pavements  supported  this  hypothesis. 

Sih's  theory  of  fracture  (10)  based  on  the  field  strength  of 
the  local  strain-energy-density  was  utilized  to  analyze  mixed 
mode  crack  propagation.  The  two  fundamental  hypotheses  of  crack 
extension  in  Sih's  theory  are: 

1.  The  crack  will  spread  in  the  direction  of  maximum 
potential  engergy  density  or  minimum  strain  energy 
density 

2.  The  critical  intensity  Scr  of  this  potential  field 
governs  the  onset  of  rapid  or  brittle  crack  propagation 


In  those  cases  where  a  fracture  is  not  a  rapid,  unstable 
process  (i.e.,  the  stress  intensity  factor  under  the  applied  load 
condition  does  not  exceed  the  critical  stress  intensity  factor, 
or  the  strain-energy-density  factor  Smin  is  less  than  the 
critical  value  Scr),  slow  stable  fatigue  crack  growth  is 
presumed.  Typically,  crack  growth  laws  relate  the  rate  of  change 
of  crack  length  to  the  stress  level  or  stress  intensity  factor 
such  as: 

dc  n  (3) 

—  -  A(AK) 
dn 
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'.V.V. 


v  v 


where 

dc/dn  *  rate  of  crack  growth 
AK  ■  stress  intensity  factor 
A,n  *  material  constents 


For  mixed  mode  fracture,  the  OSU  model  utilizes  the  crack 
growth  law  in  terms  of  the  strain-energy-density  factor  along  the 
direction  of  fracture  (Smin): 

dc  n 

—  =  B  (A  Smin)  (4) 

dn 

The  fatigue  life,  or  number  of  load  applications  to  produce  a 
crack  through  the  overlay,  is  given  by: 

cf  dc 

■  /  - 

CO  B  (ASmin) 


=  initial  starter  flaw  length 

=  crack  length  at  which  the  overlay  is 

considered  failed  (either  its  thickness  or 
the  length  at  which  the  critical  Smin  =  Scr, 
is  reached,  whichever  is  less) 

and  Scr,  B,  and  n  are  material  constants  derived  from  fatigue 
tests  on  asphaltic  concrete  beams. 

The  OSU-Fracture  Mechanics  model  is  not  a  complete  method 
for  predicting  the  occurrence  of  reflection  cracking.  An 
analytical  method  for  computing  stress  intensity  factors  and  Smin 
(such  as  a  finite  element  model)  must  be  coupled  to  a  program  to 
calculate  the  fatigue  life  in  an  incremental  fashion  using  the 
growth  law  shown  in  equation  (4).  It  is  likely  that  a  nomograph 
procedure  could  be  developed  from  this  model  similar  to  that  of 
Majidzadeh,  et  al  (11)  for  fracture  mechanics  prediction  of  load 
associated  fatigue  cracking  in  flexible  pavements.  Thus,  further 
development  of  the  OSU-Fracture  Mechanics  model  is  necessary 
before  it  could  be  implemented  by  pavement  engineers. 


Nf 

where 

co 

cf 


2.2.4  TEXAS  TRANSPORTATION  INSTITUTE  -  Fracture  Mechanics  Model 


The  TTI  model  (12)  also  uses  fracture  mechanics  crack 
propagation  theory  to  predict  cracking.  Only  Mode  I  fracture, 
and  therefore  the  K1  stress- intensity  factor,  induced  by 


Schapery's  theory  on  crack  growth  in  viscoelastic  materials  to 
develop  the  following  growth  law  (12): 


where 


and 


dc 

—  -  B 

dn  t 


2(l+l/m) 

CAK) 


[/*  wU,2h*I/n.)  dt] 


Poisson's  ratio 


maximum  tensile  stress  the  asphalt  concrete 
mixture  can  sustain 


At 

W(t) 

m 


a  dimensionless  integral  between  0  and  2 

the  period  of  the  load  cycle 

wave  shape  of  the  stress  intensity  factor 


(6) 


(7) 


slope  of  the  straight  line  portion  of  the 
tension  creep  compliance  curve  for  the  asphalt 
cement  binder 


D2  -  intercept  of  straight  line  with  log  t 
creep  compliance  curve 


0  on 


fracture  energy  density  (force  times  displacement) 
to  produce  a  unit  area  of  crack  surface. 


The  TTI  model  is  also  not  a  complete  procedure  for 
predicting  the  occurrence  of  reflection  cracking.  It  is  merely  a 
technique  for  obtaining  crack  growth  laws  without  having  to 
perform  fatigue  tests.  Fatigue  life  is  then  obtained  by 
integrating  equation  (6)  from  the  limits  of  Co  to  Cf,  similar  to 
that  of  equation  (5)  in  the  OSU  fracture  mechanics  model.  The 
limitations  discussed  for  the  OSU  fracture  mechanics  model  also 
apply  to  TTI  model. 

2.2.5  RII  -  Phenomenological  Model 

Resource  International  engineers  have  developed  a 
phenomenological  model  for  crack  prediction  in  overlaid  flexible 
pavement  structures  which  are  reinforced  by  placement  of 
engineering  fabrics  on  the  existing  surface  prior  to  overlay 
(13).  The  model  considers  only  traffic  load  stresses  in 


predicting  the  fatigue  life  of  the  overlayed  pavement;  it  has 
been  converted  into  a  computerized  design  program  called  HWYPAV. 


This  model  was  formulated  after  extensive  laboratory  testing 
had  been  conducted  which  established  the  relationship  between  the 
fatigue  life  of  reinforced  and  normal  or  unrelnforced  asphalt 
concrete  beams.  All  fatigue  tests  were  with  beams  on  an  elastic 
foundation  tested  at  a  temperature  of  70  degrees  F  (21  degrees 
C).  The  performance  factor  of  the  fabric  in  enhancing  fatigue 
life  and/or  delaying  reflective  cracking  is  called  fabric 
effectiveness  factor  (FEF),  shown  schematically  in  Figure  3: 

Nf  reinforced  (8) 

FEF  - - 

Nf  unreinforced 

i.e.,  FEF  is  simply  the  ratio  of  fatigue  lives  as  obtained  from 
the  beam  tests.  FEF  generally  ranges  from  4  to  8  depending  on 
strain  level,  placement  depth  within  the  beam,  and  fabric  type. 


The  FEF  function  is  expressed  as 

a2 

FEF  .  al  Uh)  .GEO  (9) 

where 

al  and  a2  ■  constants  depending  upon  fabric  type 

*h  *  horizontal  strain  at  the  bottom  of  the 
existing  asphalt  bound  layer. 

GEO  =  geometry  factor  which  considers  depth,  as 
shown  in  Figure  4. 

The  fat i gue  1 i f e  of  the  pavement  i n  the  HWYPAV  program  i s : 

Nf  x  Nf 'u  (FEF)  (10) 

where  Nf'u  is  a  strain  dependent  distress  function  for  asphalt 
concrete  developed  from  AASHO  Road  Test  Data 


Cracking  of  the  existing  pavement  is  accounted  for  by 
reducing  the  elastic  modulus  of  this  layer.  The  design  program 
uses  the  elastic  multilayer  program  ELSYM-5  (14)  to  calculate 
pavement  strains.  The  layers  used  in  the  strain  analyses  are 
shown  in  Figure  5.  As  can  be  seen  from  this  figure,  the  fabric 
has  not  been  included  as  a  separate  layer.  Although  the  presence 
of  fabric  has  a  significant  effect  on  the  allowable  strains,  its 
inclusion  in  layer  theory  analysis  is  of  no  practical 
significance  since  the  fabric  is  very  thin  in  comparison  with 
other  layers  in  the  system. 


The  RII  model  has  the  following  limitations: 

(a)  Since  the  model  is  phenomenological,  the  mechanics  of 
crack  propagation  and  crack  arrest  are  not  identified 

(b)  FEF  parameters  are  established  from  small  beam  tests 
and  do  not  necessarily  represent  those  of  full  slab 
condi t -  on s. 


CHAPTER  I  1  I 


MECHANISTIC  MODELS  FOR  GEOTEXTILE  SYSTEMS 


3.1  I ntroduct i on 

Field  and  laboratory  data  have  Indicated  that  fabrics 
imbedded  in  pavement  structures  enhance  the  fatigue  life  and 
retard  reflection  cracking.  To  support  these  significant 
experimental  data,  various  mechanistic  models  have  been 
postulated  which  present  a  theoretical  basis  for  the  observed 
phenomena.  The  so-called  "fabric  effectiveness"  is  attributed  to 
three  mechanisms  which  might  act  singly  or  simultaneously  to 
effectively  retard  crack  growth  in  pavements. 


3.1.1  Crack  Tip  Blunting 


To  explain  the  blunting  mechanism  of  geotextiles,  one  should 
consider  the  current  state  of  practice  in  which  bitumen-saturated 
geotextile  is  laid  on  the  top  of  an  existing  pavement  surface  or 
is  imbedded  within  an  asphaltic  layered  system. 


I n  the  case  of  rehab i 1 1 tat i on  of  an  ex i st i ng  pavement ,  as 
shown  in  Figure  5,  the  tips  of  all  existing  cracks  and 
discontinuities  in  the  old  pavement  lie  directly  beneath  the 
bi tumen- saturated  fabric.  In  the  case  of  a  new  pavement,  where 
the  geotextile  is  imbedded  in  the  lower  third  of  the  structure,  a 
propagating  fatigue  crack  tip  encounters  the  bitumen-saturated 
layer.  From  the  fracture  mechanics  viewpoint,  it  is  postulated 
that  the  crack  growth  is  a  consequence  of  the  changing  of  the 
crack  tip  profile.  Ouring  a  cyclic  deformation  of  a  pavement 
under  many  loads,  a  crack  tip  undergoes  the  phenomenon  of 
blunting  and  resharpen i ng.  The  formation  of  a  plastic  zone,  and 
its  spread  ahead  of  a  crack  tip  during  tensile  loading  cycles, 
blunts  the  propagating  crack.  Conversely,  during  unloading 
cycles,  elastic  contraction  of  the  material  surrounding  the  crack 
imposes  a  residual  stress  which  resharpens  the  crack  tip,  aiding 
its  growth.  The  rate  at  which  a  crack  may  propagate,  and  the 
path  it  follows,  depends  entirely  on  the  energy  balance  at  the 
crack  tip.  The  presence  of  a  viscous  layer  of  finite  thickness 
("bitumen-saturated  geotextile")  produces  a  large  amount  of 
plastic  deformation  which  could  blunt  the  tip  and  thus  retard 
crack  propagation.  It  is  postulated  that  bitumen- impregnated 
geotextiles,  of  finite  thickness,  could  alter  the  energy  balance 
at  the  crack  tip: 
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V  +  T  D 
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LAYER  3  EXISTING  GRANULAR  BASE 
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Figure  5.  Fabric-Reinforced  Asphalt  Overlay 
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U  “  the  work  done  by  the  external  loads 
V  =  the  elastic  stored  energy 
T  =  the  kinetic  energy,  and 

D  *  irreversible  energy,  viscous  dissipation,  plastic  work 
and  surface  energy. 

The  dot  in  the  above  quotation  indicates  the  rate  of  change  or 
differentiation  with  time.  Obviously,  the  viscous  dissipation 
rate  of  fabric-reinforced  systems  is  dependent  on  the  thickness 
of  the  tack  coat  as  well  as  the  viscoelastic  nature  of  the  fabric 
itself.  A  soft  non-woven  fabric,  with  intermediate  modulus,  when 
saturated  with  sufficient  quantity  of  viscoelastic  tack  coat, 
could  enhance  the  blunting  mechanism  as  discussed  previously. 

A  double  layer  system,  as  shown  in  Figure  6,  could  be  formed 
to  sandwich  the  viscous  material  to  form  a  viscoelastic 
constraint . 

3.1.2  Strain  Reduction 

It  has  also  been  hypothesized  that  the  addition  of 
geotextiles  reinforces  the  pavement  structure  by  increasing 
relative  stiffness  and  subsequently  reducing  the  local  strains 
and  stresses  responsible  for  fatigue  crack  propagation.  It  could 
be  argued  from  the  fracture  mechanics  point  of  view  that  a  low 
stress  intensity  factor  is  achieved  by  lowering  stresses  and 
strains  at  localized  crack  tip  regions,  which  subsequently 
reduces  the  crack  growth  rate.  Both  theoretical  and  experimental 
data  suggests  that  in  soft  ground  reinforcement  application,  the 
fabric  stiffness  may  play  a  significant  role.  However, 
unpublished  proprietary  work  recently  done  by  the  authors  has 
shown  that  for  moduli  values  in  excess  of  50,000  psi  (345  MPa) 
the  stress  intensity  factor  is  not  sensitive  to  the  stiffness  or 
modulus  of  the  bitumen-impregnated  geotextile.  Therefore ,  the 
significance  of  geotextile  stiffness,  when  imbedded  in  a  flexible 
pavement  and  saturated  with  asphalt  cement,  is  questionable.  In 
a  soft  ground  reinforcement  application,  however,  the  geotextile 
responds  as  a  tensile  membrane  with  substantial  elastic-plastic 
yielding  and  shear  distortion  occurring  in  the  soil  beneath  it. 
Theoretical  elastic-plastic  analysis  of  the  soft  ground  support 
condition  was  utilized  to  formulate  Figure  7.  This  figure  shows 
that  the  deflection  under  load  is  significantly  reduced  as  the 
fabric  modulus  increases. 


3.1.3  Buffer  Zone  Debonding 

According  to  the  Buffer  Zone  concept,  as  a  propagating  crack 
in  a  rather  stiff  medium  reaches  a  relatively  softer  zone,  the 
crack  tip  blunts  and  may  turn  its  direction  by  90  degrees, 
growing  horizontally  along  the  soft  zone.  A  crack  running  along 
the  interface  would  result  in  the  debonding  of  the 
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Fabric 


Bitumen 


Figure  6.  Double  "Sandwich"  Fabric  Layer 


asphalt -impregnated  geotextile  from  the  existing  pavement 
surface.  Recent  field  data  from  South  Africa,  using  sand  as  a 
bondbreaker  as  well  as  theoretical  analysis  of  reflection 
cracking,  confirms  that  the  debonding  mechanism  retards  crack 
growth  in  overlay  pavements,  thus  indicating  a  reduction  in  the 
stress  intensity  factor.  Factors  affecting  debonding  mechanisms 
are  fabric  thickness,  saturated  geotextile  modulus,  and  required 
modulus  ratio  between  the  asphalt  mixture  and  the  saturated 
geotextile. 


3.2  Thermal  Forces  in  Overlaid  PCC  Pavement  Structures 

It  is  generally  agreed  that  the  most  serious  reflection 
cracking  problem  occurs  in  bituminous  overlays  of  PCC  pavements. 
It  is  also  generally  agreed  that  the  primary  cause  of  the 
cracking  is  thermally  induced  movement  of  the  concrete  slab.  The 
literature  review  of  theoretical  models  has  recently  disclosed 
much  disparity  regarding  the  stress  state  in  the  overlay  under 
thermal  loading. 

Thermal  stresses  result  from  both  seasonal  and  daily  changes 
in  slab  temperature.  The  thermal  loading  can  be  represented  by 
the  superposition  of  two  different  thermal  conditions: 

A.  Uniform  change  (  T)  in  slab  temperature.  This  condition 
represents  seasonal  changes  in  average  slab  temperature 
which  can  occur  over  long  time  periods. 

B.  "Pure"  curling.  This  condition  represents  the  daily  or 
short  time  period  temperature  variation  within  the  slab. 
For  pure  curling  the  average  slab  temperature  has 

not  changed;  however,  the  top  of  the  slab  is  colder  than 
the  bottom  of  the  slab  with  the  temperature  assumed 
linearly  related  to  slab  depth.  The  curling  gradient 
(CG)  is  given  in  degrees  F/in  (degrees  C/mm)  of  slab 
depth.  Figure  8  shows  the  representation  of  thermal 
loading  using  these  two  definitions.  The  reference 
temperature  (TR)  is  the  "zero-stress"  temperature  for  the 
overlay;  slab  temperatures  below  TR  will  transfer 
tensile  stresses  to  the  overlay.  Figure  9  shows  the 
expected  monthly  average  slab  temperature  and  curling 
gradient  for  overlaid  concrete  slabs  in  Ohio.  This 
figure  is  based  upon  computer  prediction  of  pavements  in 
the  central  Ohio  area.  Slab  thickness  varied  from  8  to 
10  inches  (203  to  254mm)  and  asphalt  overlay  thickness 
from  2.5  to  5.0  inches  (64  to  127mm)  for  the  pavements 
used  in  the  Ohio  study.  Figure  9  provides  an  estimate  of 
the  thermal  load  magnitudes.  Expected  curling  gradients 
(CG)  varies  from  .5  degree  F/in  (.011  degree  C/mm)  in  the 
spring  and  fall  to  about  1  degree  F/in  (.022  degree  C/ram) 
in  the  winter  months.  Mean  slab  temperature  changes  by 
about  40  degrees  F  (22  degrees  C)  from  summer  to  winter, 
dropping  at  a  rate  of  about  8  degrees  F  (4.4  degrees  C) 
per  month  during  the  fall. 


o 

(1) 

Li 

3 

<U  3  ±J 

- — * 

■U 

>  4J  C 

• 

C  J2 

fl  4) 

c 

ru  m 

Li 

Cu 

•U  Li  rl 

•H 

o;  rH 

O 

0 

•-»  4»  T) 

\ 

£  V) 

a 

U)  a  i0 

Cu 

e 

0  e  Li 

0 

tv 

0i  DO 

e-> 

E- 

Expected  Seasonal  Pavement  Temperature  and  Curling  Temperature  Gradient 


The  2-D  finite  element  analysis  of  the  full  scale  pavement 

shown  in  Figure  10  was  conducted  using  the  SAP  IV  program.  As 
shown  in  Figure  10,  only  asphalt  overlay  modulus  and  overlay 
thickness  were  varied  in  these  analyses.  Slab  length  was  20  feet 
(6.1m)  in  all  cases.  Two  separate  thermal  loading  conditions 
were  analyzed:  (1)  a  uniform  reduction  (  T)  of  30  degrees  F 
(16.7  degrees  C)  in  slab  temperature,  and  (2)  "pure"  curling  with 
gradient  (CG)  of  .5  degree  F/in  (.011  degree  C/mm).  For  constant 
overlay  thickness  and  modulus,  overlay  stresses  and  joint  opening 
were  found  to  be  linearly  related  to  T  and  CG.  Both  full 
friction  (no  slip)  and  no  friction  (slip)  between  PCC  slab  and 
aggregate  base  were  investigated.  Full  friction  reduced  overlay 
maximum  stress  by  less  than  6%  for  uniform  temperature  change  and 
4%  for  curling  as  compared  to  the  no  friction  condition.  Full 
bond  between  asphalt  overlay  and  PCC  slabs  was  assumed  in  all 
cases . 


Figures  11  through  14  show  computed  stresses  in  the 
overlay  at  the  center  of  the  joint  as  a  function  of  depth  (Z) . 
In  all  cases  for  both  curling  load  and  uniform  temperature 
change,  maximum  stress  occurs  at  the  bottom  of  the  overlay.  The 
stress  distributions  are  very  similar  for  the  uniform  temperature 
change  and  curling  loading  conditions.  The  similarities  occur 
throughout  the  range  of  overlay  thickness  and  overlay  moduli 
investigated.  Figures  15a  and  b  show  the  computed  shear  stress 
at  the  overlay/slab  interface.  Again  the  stress  distributions 
are  similar  for  the  two  loading  conditions.  These  shear  stresses 
would  have  to  exceed  the  tack  coat  bonding  stress  to  cause 
slippage  between  the  two  layers.  The  maximum  shear  stresses  are 
below  the  bonding  stresses  given  by  ARE  (8) .  Figure  16  shows  the 
effect  of  breaking  the  bond  (either  by  slippage  or  introduction 
of  a  bond  breaker)  upon  maximum  overlay  stress.  A  dramatic 
reduction  in  stress  is  predicted  for  bond  breaker  lengths  as 
short  as  1  in  (25mm) . 


The  high  sensitivity  of  overlay  stress  to  overlay  stiffness 
is  clearly  demonstrated.  As  noted  earlier,  ARE  (8)  suggests  that 
the  creep  modulus,  Ec,  be  utilized  for  stress  calculations. 
However,  Ec  is  both  temperature  and  time  of  loading  dependent. 
Figure  17  presents  this  dependency  for  a  typical  dense-graded 
asphalt  concrete  with  Ec  calculated  by  the  Heukelom  and  Klomp 
(15)  (or  Shell)  procedure.  This  procedure  determines  the 
compressive  creep  modulus.  The  tensile  creep  modulus  is  actually 
needed  for  reflection  cracking  analysis.  However,  no  procedure 
for  predicting  tensile  creep  modulus  has  been  published.  If 
creep  modulus  curves  in  tension  are  similar  to  those  of  Figure 
17,  then  the  implications  to  thermal  reflection  cracking  analysis 
and  modeling  are  very  significant.  An  incremental  analysis  which 
utilizes  the  loading  time  and  temperature  dependent  creep  modulus 


Figure  10.  Schematic  2-D  Finite  Element  Model 
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Figure  15a.  Shear  Stress  Overlay  Due  to  Uniform  Change  in  Slab  Temperature 
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Figure  15b.  Shear  Stress  Distributions  Overlay  Due  to  "Pure"  Curling 
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Figure  16.  Effect  of  bond  breaker  length  or  stresses 
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Figure  17.  CREEP  MODULUS  VS.  TEMPERATURE 


would  be  necessary  to  calculate  overlay  stress  and  joint  opening. 
Seasonal  changes  occur  over  long  periods  of  time  (time  required 
to  drop  from  TR  to  T  amount}  while  curling  can  occur  over 
relatively  short  periods  (less  than  1/2  day)  and  at  all 
temperatures.  The  fact  that  curling  occurs  over  shorter  loading 
times  than  seasonal  uniform  temperature  change  means  that  a 
higher  Ec  should  be  used  for  curling  load  than  for  uniform 
temperature  change  stress  calculations.  The  higher  Ec  will 
result  in  higher  overlay  stresses  which  should  be  considered  when 
comparing  seasonal  and  curling  Induced  loading  conditions. 

Since  both  curling  and  seasonal  change  produce  joint 
openings  which  induce  similar  overlay  stress  distributions,  it  is 
believed  that  it  is  possible  to  equate  the  two  loading 
conditions.  Figure  18  shows  that  overlay  stress  is  linearly 
related  to  displacement  at  the  top  of  the  joint.  The  slope  of 
the  stress- joint  displacement  line  is  a  function  of  overlay 
modulus.  The  overlay  stresses  are  Independent  of  the  type  of 
thermal  load.  At  constant  overlay  modulus,  a  horizontal 
displacement  of  x  will  produce  the  same  overlay  stress  regardless 
of  whether  this  displacement  was  produced  by  slab  curling  or 
uniform  temperature  change.  This  is  further  evidence  that 
curling  can  be  equated  to  seasonal  temperature  change. 

Figure  19  shows  the  relationship  between  overlay  modulus  and 
predicted  joint  displacement  for  the  finite  element  model.  For 
Ec  below  1500  psi  (10.3  MPa),  movement  nearly  equals  that  for 
free  unrestrained  thermal  movement  (  TL/2).  From  Figure  17  the 
expected  modulus  is  below  1500  psi  (10.3  MPa)  for  long  loading 
times  (greater  than  1  month)  for  temperature  greater  than  25 
degrees  F(-4  degrees  C) .  Subsequent  design  of  fabric  reinforced 
overlay  should  consider  that  the  modulus  of  the  overlay  will  be 
less  than  1500  psi  (10.3  MPa).  The  joint  opening  can  then  be 
calculated  using  the  a&T  L/2  expression.  This  will  simulate 
joint  openings  which  occur  in  real  pavements  at  temperatures 
above  about  25  degrees  F  (-4  degree  C). 
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Figure  19.  Thermal  Movement  of  slab  as  a  function  of  asphalt  creep 
modulus. 
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CHAPTER  IV.  DESIGN  PROGRAM 


4.1  Scope  and  Outline  of  Research 


The  primary  objective  of  this  study  is  to  develop  a 
procedure  for  the  design  of  flexible  overlays  for  existing  rigid 
pavements  based  upon  mechanistic  concepts  of  the  "Overlay 
Pavement  Foundation"  (OPF)  structure,  including  the  use  of 
engineering  fabrics  as  a  reflective  crack  arrest  phenomenon. 


The  program  can  also  be  used  to  design  rigid  overlays  of 
rigid  or  flexible  pavements,  as  well  as  flexible  overlays  of 
flexible  pavements,  if  appropriate  material  properties  are  input. 


A  three-dimensional  finite  element  approximation  would 
probably  be  most  rational  representation  of  the  OPF  system. 
However,  the  amount  of  linear  equations  involved  in  the  solution 
of  any  non-academic  problem  is  so  large  that  it  makes  this 
approx i mat i on  economically  unreasonable  and  very  sensitive  to  a 
numerical  error.  Therefore,  in  this  study  a  two-dimensional 
model  is  used  to  determine  the  stress  patterns  in  the  overlay 
structure  and  a  three-dimensional  model  (RISC*) computer  program 
(16))  can  be  used  to  estimate  the  existing  pavement  life. 


4.1.1  Model  Description  and  Program  Capacity 


For  finite  element  modeling  the  OPF  structure  is 
into  the  following  structural  elements  (see  Figure  20): 
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Each  of  these  structural  elements  has  its  own  finite  element 
type,  mesh  and  materials  properties,  which  are  described  in  the 
following  sect i ons . 
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4.1.2  Overlay  Structure 


For  finite  element  approximation  of  the  overlay  structure  a 
two-  dimensional  1 soparametr Ic  element  with  incompatible  mode 
(Wilson  element)  are  used  (see  Figure  21).  The  local  and  global 
coordinates  are  related  by  the  transformat  Ion 
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1*1,  2,  3,  4,  are  the  global  displacement  components 
the  global  x  and  y  directions  respectively. 

To  account  for  the  errors  in  bending  deformation,  Wilson 
(17)  suggested  that  the  following  form  of  displacement 
approximation  be  used  if  higher  accuracy  is  desired. 


Elastic 

Foundation 


V 

Existing 

Slabs 


Figure  20.  Structural  dominates  of  the  pavement 


u  ■  \f/  u 


(19. a) 


T 

v  *  \ff  v 


(19. b) 


where 


[V  u2*  V  V  V  a23 


(20. a) 


v  =  tv,,  v2,  v3,  v4,  a3,  a4] 


(20. b) 


^ ( I-s)  (l-tA 
(1+8)  (1-t) 

(1+S)  (1+t) 

(1-S)  (1+t) 


+  *  1/4 < 


•I 


4  (1-s  ) 


* 

V 


(21) 


^  4  (1-t  )  ) 

/ 

(1-S  ) 


(22) 


(1-t  ) 


Here  the  a's  are  arbitary  constants  representing  internal  degrees 
of  freedom  within  the  element.  Hence  the  displacement  matrix  can 
be  written  as 
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where  the  matrices  <#>  and  X  are  given  by  Equations  (16)  and  (22). 

4.1.3  Strain-Displacement  Equations 

For  a  two-dimensional  analysis,  the  strain-displacement 
equations  are  given  by: 
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In  view  of  Equations  (19),  (20)  and  (21),  Equation  (28)  becomes 


where  ^  fs  the  transpose  of  'J'  and  'J',*  are  partial 

derivatives  of  ^  with  respect  to  x  and  y,  ,y 
respectively.  With  the  chain  rule  it  can  be  verified 
that 
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where  the  Jacobian  determinant  j|  is  given  by 
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for  the  plane  strain  condition,  and  for  the  plane  stress 
cond i t i on 
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=  modulus  of  elasticity  in  the  t  direction 


(40) 


=  Poisson  coefficient  which  characterizes  the 
decrease  In  the  j  direction  for  tension  in 
i  d  I  rect  i  on 

=  shear  modulus  in  the  x-y  plane. 


Here  it  is  assumed  that  the  principal  axes  for  the  orthotropic 
material  coincide  with  the  global  axes.  If  this  is  not  the  case, 
the  following  transformation  has  to  be  made.  Let  be  the  angle 
of  rotation  from  the  global  x-y  axes  to  the  material  principal 
x'-y'  axes.  Let  C  be  the  compliance  matrix  with  respect  to  the 
material  principal  axes  and  C  be  the  corresponding  compliance 
matrix  with  respect  to  the  global  x-y  axes;  the  third  material 
principal  axis  (normal  to  the  x'-y'  plane)  being  coincided  with 
the  global  z  axis.  Then  the  two  matrices  are  related  by; 
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Having  obtained  C  in  the  global  coordinate  system,  the  matrix  E 
can  be  computed  as  the  inverse  of  the  C  matrix  and  substituted  in 
Equation  (38).  The  Integrals  in  (38)  are  carried  out  by  the 
direct  application  of  one-dimensional  numerical  integration  using 
the  two-point  Legendre  integration  formula.  For  example. 


Figure  22.  Principal  Axes  of  Material  Properties  and  Global  Axes 
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where  the  Integration  points  sj,  ti,  and  the  corresponding  weight 
functions  are  given  by 

s  =  t  m  +  0.57735026918963 

1  1 

s  =  t  =  -  0.57735026918963 

2  2 

(44) 

W  =  W  =1.0 
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4.1.5  Element  Equivalent  Nodal  Load  Vector 


(a)  Due  to  the  Body  Force 


Since  the  body  force  considered  here  is  due  to  gravity,  it 
is  assumed  to  be  constant  throughout  the  element.  Let  fx  and  fy 
be  the  components  of  gravity  force  in  the  global  x  and  y 
directions  respectively,  i.e.. 


(45) 
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(b)  Due  to  Temperature 

The  temperature  variation  within  the  element  Is  assumed  to 
be  related  to  the  temperatures  at  the  nodes  of  the  element  by  the 
Interpolation  functions^,  t.e.. 


/  T 


where 


Ft  ,T  ,T  ,T  1 

l  1  2  3  4j 


temperature  at  node  1 


(48) 

(49) 
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<t>  is  given  by  Equation  (16).  The  temperature  rise  at  any  point 
within  the  element  Is  assumed  to  be 


AT 


T  -  T. 


(50)* 


T  «  temperature  rise  at  any  point  (s,t) 

r  *  reference  temperature  for  the  temperature 
r  stress  free  state  within  the  element 


For  an  orthotropic  material ,  it  can  be  shown  that  the 
initial  strains  due  to  temperature  change  AT  is  given  by; 


(51) 


for  plane  strain  condition,  and  for  plane  stress  condition,  by* 


(52) 


where 
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If  the  principal  axes  of  the  material  do  not  coincide  with  the 
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global  axes,  the  strain  transformat  Ion  has  to  be  performed,  using 
the  transformat Ion  matrix  T77  given  in  Equation  (41),  in  the  same 
manner  as  the  compl iance  matrix  transformat ion  (40).  The 
equivalent  nodal  load  vectors  due  to  temperetqre  change  are: 

P  s  b  f  f  |j|  B  E  e  ds  dt  (53) 
~q  l\  J-l  '  q  ~0 


P  a  b  f  I  |j|  B  E  €  ds  dt  (54) 

~r  -L|  r  ^0 

The  integrations  are  carried  out,  as  before,  by  the  two-point 
L  egendre  i nt egrat i on  f ormu 1  a . 

(c)  Due  to  Boundary  Surface  Tractions 

For  an  element  which  has  a  boundary  face  subjected  to  a 
distributed  external  load,  the  equivalent  nodal  loads  due  to  this 
surface  load  have  to  be  considered.  For  illustration  purposes, 
it  is  assumed  that  the  boundary  1-2  of  the  element  is  subjected 
to  the  prescribed  linearly  distributed  load  per  unit  area  as 
shown  below. 


Figure  23.  An  Element  Subjected  to  Surface  Pressure  on  1-2  Side 


Let  pxl,  pyl  be  the  components  of  the  load  intensity  pi  at 
node  1  and  px2,  py2  be  the  components  of  the  load  intensity  p2  at 
node  2.  Since  along  1-2  edge  the  local  t  ordinate  is  constant  («= 
-l),  hence  it  can  be  shown  that  the  components  px  and  py  of  the 
distributed  load  along  this  surface  are  described  as: 
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Along  the  1-2  side,  the  matrix  of  displacement  interpolation 
f unct i ons $ q  and$r  are  given  by 
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t=  -1 

The  differential  surface  dS  along  the  side  1-2  is  given  by: 


dS  =  b 


i?  ♦  l*\ 2 


ds  =  b l  ds 
12 


where  l  is  the  length  of  side  1-2  and  b  is  the  element 
12 

constant  thickness. 

The  reduced  equ i va 1 ent  noda 1  1 oad  vector  i s 

T  -1 

P  .  p  -  K  K  £  (59) 

~  ~q  rq  rr  r 

In  the  case  when  the  incompatible  modes  of  displacements  are  not 
used, 

P  »  P  (60) 


If,  in  addition,  a  set  of  concentrated  loads  are  applied  at  the 
nodes,  the  final  element  load  vector  will  be  of  the  form: 


q  rq  rr  r  w 

where  £  N  is  the  vector  of  concentrated  load  at  the  nodes 
(or  the  static  equivalent  of  a  simple  beam 
reactions  of  a  distributed  load  between  the 
nodes) . 


4.1.6  Element  Stress  Output 


The  element  output  stress  vector  is  given  by: 

tr  ™  E  (c  -  c  ) 

~  ~0 


[°rx,<ry,rxy] 


«,  can  be  expressed  as : 


-1  -1 

«  -  (B  -  K  K  )  q  +  B  K  P  (64) 
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Hence,  Equation  (62)  can  be  written  as: 

-1  -1 

<r  **  E  (B  -  K  K  )  q  +  E(B  K  P  -e  )  (65) 
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Once  the  solution  for  the  compatible  nodal  displacement  vector  q 
is  obtained,  the  element  stresses  at  some  specified  points  in  the 
element  can  be  computed  from  Equation  (65).  Obviously,  if  the 
incompatible  modes  of  displacement  are  not  used,  Equation  (65) 
reduces  simply  to 

cr  ■  EB  q  -  E  e  (66) 

^  'v  ^  r\ 

q  0 

Here, ~0  is  the  initial  strain  vector  which,  if  due  to 
temperature  changes,  is  given  by  either  Equation  (51)  or  (52). 

4.1.7  Material  Properties 

The  overlay  structure  is  assumed  to  be  a  homogeneous  and 
isotropic  matrix;  therefore,  the  elasticity  matrix  for  plane 
stress  condition  is  given  by: 
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L  symmetric  j 
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and  for  plane  strain  condition  is  given  by: 
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4.2  Fabric 


The  fabric  in  this  computer  model  has  been  represented  as  a 
N- layered  stripe  (N  =  0,1,. ..,4).  Each  layer  could  have 
individual  orthotropic  material  properties,  thicknesses  apd  .fabric 
lengths  (see  Figure  24).  -v 


With  this  model  it  is  possible: 


(1)  To  investigate  the  Fabric  Effectiveness  Factor  (FEF)  of 
the  composite  fabric  structure. 


(2)  To  model  bond  conditions  between  fabric  and  overlay  or 
f abr i c  and  ex i st i ng  s 1 ab . 


(3)  To  establish  a  necessary  fabric  length  in  the  case  of 
horizontal  temperature  change  when  fabric  must  act  as  a 
bond  breaker  to  prevent  thermal  cracking  of  the  overlay. 


Finite  elements  used  to  approximate  the  fabric  absolutely 
are  identical  to  "overlay  elements"  with  the  exceptions  that 
Equations  (39)  and  (40)  are  applicable. 


4.3  Existing  Slabs 

Existing  slabs  generally  have  the  same  finite  element  type 
as  the  overlay  structure  described  above  except  for  the  elements 
used  for  modeling  cracks.  These  elements  have  a  very  small 
elastic  modulus  In  x,y,z  directions  and  Poissons'  ratio  vx  =  vy  * 
•'z  »  .47.  In  other  words  they  represent  a  crack (s)  with  no  shear 
resi stance. 


4.4  Slab-Foundat Ion  Contact  Zone 


The  correct  finite  element  model  of  this  part  of  the  OPF  is 
extremely  important  In  the  case  of  a  temperature  load  or  a 
vertical  mechanical  load  which  can  produce  either  a  partial 
contact  of  the  slab  or  a  large  horizontal  df splacement.  Partial 
contact  due  to  a  void  formation  will  be  discussed  later.  Special 
"Joint"  elements  are  used  in  this  model  to  achieve  more  realistic 
finite  element  approximation  of  the  contact  zone.. 

The  Joint  element  was  first  introduced  by  Goodman  (18)  as  a 
finite  element  model  representing  joints  and  seams  in  the 
analysis  of  rock  mass  structures.  The  configuration  of  the  joint 
element  is  as  shown  in  Figure  25.  It  consists  of  a  pair  of 
straight  lines  with  four  nodal  points.  The  element  has  zero 
width,  a  unit  length,  and  a  thickness  b  in  the  direction  normal 
to  the  plane  of  the  element  (x-y  plane).  The  nodal  point  pairs 
(1,  4)  and  (2,3)  initially  have  identical  coordinates. 

The  joint  element  is  assumed  to  have  essentially  no 
resistance  to  a  net  tension  in  the  direction  normal  to  the  sides 
1-2  and  3-4,  i.e.,  in  the  y'  direction.  It  offers  high 
resistance  to  compression  in  the  y'  direction  and  may  deform 
somewhat  under  normal  pressures,  particularly  if  there  are 
crushable  irregularities  or  compressible  filling  materials.  It 
also  offers  shear  resistance  in  the  x'  direction  under  normal 
pressures.  However,  at  no  point  within  the  element  can  the  shear 
stresses  exceed  shear  strength  of  the  joint. 

Linear  displacements  are  assumed  for  the  element.  Let  v' 
and  u'  be  the  displacements  in  the  normal  and  the  tangential 
directions  of  the  element,  respectively.  Then  the  d i sp 1 acement s 
along  the  top  of  the  element  can  be  expressed  in  terms  of  the 
local  nodal  displacements  as: 
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Similarly,  along  the  bottom  of  the  element 
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Next  let  Wn  and  Ws  be  the  relative  displacements  between  the  top 
and  bottom  of  the  element  in  the  normal  apd  tangential 
directions.  Hence,  in  view  of  Equations  (69)  and  (71),  the 
relative  displacements  are  given  by: 

w  »  Bq1  (72) 
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The  local  nodal  displacements  are  related  to  the  global  nodal 
displacements  by  the  transformation 

q'  **  T  q  (76) 
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Fig.  25  A  Joint  Element  in  bocal  and  Global  Coordinates 
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The  element  has  a  normal  stiffness  Kn  and  a  shear  stiffness 
Ks  (both  in  F/l_3  units).  The  stresses  at  any  point  within  the 
element  are  proportional  to  the  relative  displacements  between 
the  top  and  bottom  of  the  element,  i.e.. 
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The  element  stiffness  matrix  in  the  global  coordinate  system 
is  g i ven  by 


k  a  b  T  (/  B  k  B  dx')T 


with  the  displacement  Interpolation  matrix  B  given  by  Equation 
(75)  and  the  coordinate  transformat  ion  matrix  Tc  giver)  by 
Equation  (78).  The  axial  stiffness  k  in  Equation  (83)  is 
replaced  in  this  case  by: 


Substituting  Equations  (75)  and  (84)  Into  Equation  (83)  and 
integrating  over  the  element  leads  to 


K  =  T  K'  T 


where  K,  is  th^  global  element  stiffness  matrix,  and 
K'  is  the  local  element  stiffness  matrix  given  by: 
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The  stress  output  can  be  computed  from  Equation  (88)  which,  in 
view  of  Equation  (84)  and  (85)  can  be  written  as: 

o-  =  k  T  q  (89) 


where  Z  -  [<r  •  T]  (90) 

Thf  correspond i ng  equivalent  nodal  load  vector  is  calculated 


where 


T  [*P  ,P  ,P  ,P  ,P  ,P  ,P  ,P  ] 

P  =  |xxxxyyyy|  (92) 
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PXj  and  pyj  are,  respectively,  the  equivalent  nodal  load 
components  in  the  x  and  y  directions  at  node  i. 

The  resulting  shear  stresses  must  be  compared  with  the  shear 
strength  of  the  element.  The  shear  strength  at  any  point  within 
the  element  is  assumed  to  be  of  the  forms 

Ts  ~  c  +  or  tan  $  (93) 

where  Ts  is  the  shear  strength,  c  is  the  cohesion  and  <f>  is  the 
angle  of  friction  of  the  material,  <r  is  the  normal  compressive 
stress  at  that  particular  point.  The  shear  stresses  at  the  nodal 
points  are  checked  with  the  corresponding  shear  strength.  If  at 
any  node  the  stress  exceeds  the  shear  strength,  it  is  reduced  to 
the  shear  strength  value.  The  excess  shear  stress  is  used  to 
compute  the  correspond i ng  equivalent  nodal  loads  and  are  applied 
back  as  a  new  load  vector  in  the  iteration  procedure. 

The  element  normal  stresses  are  also  checked  for  tension. 
If  any  nodal  normal  stress  is  tensile,  the  normal  and  shear 
stresses  at  that  particular  node  are  reduced  to  zero  and  the 
corresponding  nodal  loads  are  applied  back  as  a  new  set  of  load 
vectors  in  the  iteration  procedure.  The  iteration  process  is 
repeated  until  all  of  the  nodal  normal  stresses  are  essentia.lly 
compressive  and  the  nodal  shear  stresses  do  not  exceed  the  nodal 
shear  strength.  The  iteration  scheme  used  is  based  upon  the 
initial  stress  concept  proposed  by  Zienkiewicz  (19). 


4.5  Elastic  Foundation 

A  complete  review  of  the  State-of-tiie  Art  of  existing  models 
of  the  pavement  foundations  couple  is  presented  in  Chapter  II  of 
Reference  [16];  Chapter  III  of  the  same  reference  discusses  the 
model  used  in  the  RISC  computer  program  development. 

An  identical  procedure  has-been  incorporated  in  this  study, 
with  corrections  for  a  2-dimensional  space,  i.e.,  all 
integrations  by  area  in  3-dimensional  case  have  been  replaced  by 
integration  by  line.  In  addition  to  elastic  foundation  support, 
the  EFRON  computer  program  provides,  as  an  option,  Winkler  or 
rigid  type  of  support.  This  option  is  provided  so  that  the  user 
can  compare  the  results  from  the  two  different  methods.  It 


should  be  emphasized  that  Winkler  foundation  Is  not  a  reasonable 
representation  of  a  real  foundation  and  should  not  be  used  In 
general . 

An  Initial  void  could  be  assumed  under  existing  slab,  as 
shown  in  Figure  26.  In  a  finite  element  model,  a  void  Is 
reflected  by  replacing  with  zero  all  diagonal  and  off-diagonal 
terms  of  the  fully  populated  stiffness  matrix  of  the  foundation 
which  correspond  to  the  nodal  points  out  of  contact,  eg.  (94). 


where  K.  *  stiffness  matrix  of  the  foundation 

N  »  number  of  nodal  points  along  the  slab(s) 

4.6  Joint  Structure 

Doweled  or  plain  joints  can  be  presented  in  finite  element 
model,  as  shown  in  Figure  27. 

Joint  space  Is  represented  in  the  model  by  a  set  of 
rectangular  elements  with  very  low  elastic  moduli  in  all 
directions.  The  dowel  bar,  if  it  exists,  is  approximated  by  one 
of  the  three  types  of  elements  indicated  below: 

1.  Type  1  elements  correspond  to  an  actual  body  of  the  dowel 
bar  and  is  represented  by  element  type  described  in 
Section  2.1  with  elastic  modulus  Ee  computed  from: 


where 


h 

e 


D 


D.D. 


E  ■  Elastic  modulus  of  the  dowel  bar 

D.B. 

D  *  Dowel  bar  diameter 

D.B. 

h  *  Elements  height 

e 

b  *  Elements  width 

e 


2.  Type  2  and  3  elements  actually  model  the  horizontal 
degrees  of  freedom  of  the  dowel  bar.  These  elements  have  a  very 
low  horizontal  elastic  modulus  E  .  In  addition  to  that>  the  type  2 
element  has  a  very  low  shear  modulus  G. 
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Figure  27.  Finite  element  representation  of  plain  and  doweled  joints. 


V 


The  major  reason  for  such  a  complicated  model  of  the  dowel 
bar  Is  to  transfer  bending  moment  from  slab  to  slab.  This  would 
not  be  possible  In  any  other  finite  element  model  using  bar  or 
beam  types  of  elements  because  of  the  existence  of  pnly  two 
degrees  of  freedom  In  each  nodal  point  of  the  2-dlmenslonal 
element. 

To  Incorporate  dowel  bar  looseness,  an  effective  (reduced) 
dowel  bar  diameter  Is  used,  which  can  be  calculated  from 

I 

.7 

D'  =  D  [ 1-10.5L-12.58L  1og(E*/10Q0) ]  (97) 

D.B. 

where 


0*  =  Effective  dowel  bar  diameter 
b  =  Actual  dowel  bar  diameter 
O.B. 

L  =  Amount  of  looseness  (inches) 
Es  *  Subgrade  modulus  (PS I) 


This  equation  is  p  result  of  the  parametric  study  described 
I  n  Chapter  III  of  Ref erence  (16). 


4.7  Cracks 

A  finite  element  model  of  the  crack (s)  is  absolutely 
identical  to  a  plain  joint  model  described  above.  Up  to  ten 
different  cracks  could  be  included  in  a  single  problem.  The 
crack  pattern  is  not  necessary  nor  does  it  have  to  be  identical 
for  both  s 1 abs  I f  vert  leal  1 oad  I s  present .  However  for  a 
horizontal  temperature  load  an  identical  crack  pattern  must  be 
used  if  cracks  are  present. 


CHAPTER  V.  LABORATORY  TESTING 


5.1.  Types  of  Testing 

This  study  is  concerned  with  the  evaluation  of  the  cracking 
resistance  of  the  fabric  reinforced  asphalt  concrete  overlays  of 
rigid  (Portland  cement  concrete)  pavement.  In  order  to  evaluate 
the  cracking  resistance  of  the  fabric  reinforced  system,  tests 
were  conducted  under  the  f o 1  lowing  cond i t i ons  at  two 
temperatures : 

1.  Fatigue  resistance  under  load  of  asphalt  overlays  over 
concrete  bases  for  both  control  and  fabric  reinforced 
conditions.  Control  conditions  are  conventional  overlays 
without  fabric  reinforcement. 

2.  Horizontal  and  vertical  temperature  induced  stresses: 
test  condition  for  both  control  and  fabric  reinforced 
overlays. 

The  beam  testing  simulates  full  scale  flexible  pavement 
overlay  behavior  under  aircraft  load  as  experienced  in  runways 
and  taxi ways.  Figure  28  illustrates  the  concept  for  an  overlay 
condition.  Fabric  re i nforcement  involves  placement  of  an 
engineering  fabric  or  geotextile  underneath  the  asphalt  overlay. 
The  purpose  of  fabric  re i nforcement  is  to  enhance  the  fatigue 
crack  resistance  of  the  bound  layer.  Simulation  of  reflective 
cracking  of  rigid  pavement  overlays  requires  modeling  of  both 
thermal  and  traffic  loading  conditions.  Thermal  stresses  result 
from  both  seasonal  and  daily  changes  in  slab  temperature.  The 
thermal  loading  can  be  represented  by  the  superposition  of  two 
different  thermal  conditions  as  was  discussed  in  Section  3.2. 

a.  Uniform  change  in  slab  temperature:  this  condition 

represents  the  seasonal  changes  in  average  slab 
temperature  which  occur  over  long  periods  of  time. 

b.  Pure  curling  condition:  daily,  or  short  time  period, 
temperature  variation  within  the  slab.  For  pure  curling 
the  average  slab  temperature  remains  unchanged  but  the 
upper  surface  of  the  slab  is  at  a  different  temperature 
from  the  lower  surface,  with  the  temperature  assumed  to 
be  linearly  related  to  the  slab  depth. 


Manufacturing  full  scale  pavement  models  and  subjecting  them 
to  actual  thermal  loads  experienced  in  the  field  is  neither 
economically  feasible  nor  possible  given  the  time  constraints  of 
this  study.  Therefore  model  pavements  with  external  forces 
applied  to  produce  joint  movements  equal  to  those  of  full  scale 
pavements  under  field  thermal  loading  were  designed.  Uniform 


seasonal  reduction  in  slab  temperature  was  simulated  by  applying 
horizontal  tensile  forces  to  the  PCC  slab  pavements.  Traffic 
forces  were  simulated  by  applying  dynamic  vertical  loads  to  the 
model  which  was  a  beam  supported  on  an  elastic  foundation. 


5.2  Fabric  Selection 

Three  types  of  engineering  fabric  membranes  were  selected 
for  this  study  and  a  fourth  type  of  membrane  is  currently  being 
evaluated.  The  test  results  of  this  material  will  be  included  in 
the  final  report  at  the  end  of  all  work.  The  three  fabrics  were 
selected  to  represent  typical  fabrics  with  low,  medium,  and  high 
tensile  moduli.  The  pertinent  fabric  properties  are  listed  in 
Table  1. 

5.3  Mix  Design 

A  brief  report  on  the  mix  design  of  concrete  used  for  rigid 
base,  and  mix  design  used  for  asphalt  overlays  are  detailed  in 
the  following  sections: 

5.3.1  Cement  Concrete  Beams 

Two  different  sizes  of  aggregates  from  the  local  American 
Aggregate  Company  were  utilized  for  the  manufacture  of  cement 
concrete  beam  specimens.  The  aggregates  used  were  If 57  crushed 
gravel  as  coarse  aggregate  and  natural  sand  as  fine  aggregate. 
The  aggregate  gradation  limits  were  in  accordance  with  FAA 
specifications  (401),  actual  gradation  chosen  for  preparing  the 
mix  design  is  shown  in  Figure  29.  The  quantities  of  the 
constituents  of  concrete,  namely  coarse  and  fine  aggregate, 
cement  and  water,  were  chosen  as  per  FAA  specifications  for  Class 
C  concrete  after  making  necessary  adjustments  for  water  absorbed 
by  the  dry  aggregate,  and  for  the  net  water  requirement  of 
concrete  due  to  presence  of  water  in  the  polymer  additive.  The 
polymer  additive  .  used  was  supplied  by  Dow  Chemical  Company  and 
the  quantity  used  was  10  gallons  of  additive  per  cubic  yard  of 
concrete.  The  calculated  mixture  preparations  were  checked  by 
means  of  trial  batches.  The  quantities  by  weight  of  the  several 
constituents  of  regular  concrete  are  indicated  in  Table  2.  The 
slump  selected  was  three  inches,  entrained  air  4.5%  and  water 
cement  ratio  of  0.45.  The  procedures  delineated  in  ASTM  C-192 
were  strictly  followed  for  making  and  curing  beam  specimens  of 
concrete  in  the  laboratory.  Molds  used  for  casting  beam 
specimens  were  made  out  of  sawed  wood.  Molds  were  water-tight 
during  use,  as  judged  by  their  ability  to  hold  water.  Beam  molds 
were  rectangular  in  shape  and  of  the  dimensions  required  to 
produce  the  stipulated  specimens.  Curing  of  concrete  specimens 
was  performed  as  per  ASTM  C-192.  After  the  required  period  of 
curing,  the  bleam  specimens  were  sawed  with  a  diamond  blade  saw 
to  get  specimens  of  the  following  sizes: 

(1)  1.5  x  3  x  12  in  (38  x  76  x  305  mm) 


TABLE  1 


FABRIC  PROPERTIES 


PROPERTY  FABRIC  B  FABRIC  C  FABRIC  D  TEST 

METHOD 


Weight  (02) 

Modulus  (psi) 

Tensile  Strength  (lbs) 
Elongation  (percent) 
Burst  Strength  (psi) 
Asphalt  Retention  (gsy) 
Shrinkage 


3.8 

5. A 

3 

NA 

NA 

850 

ASTM  D  1682 

90 

90 

90 

ASTM  D  1682 

55 

55 

55 

ASTM  D  1682 

NA 

230 

NA 

Mullen  Burst  Test 

0.20 

0.18 

0.07 

Texas  DOT-3099 

NA 

<2 

.3 

Texas  DOT-3099 

TABLE  2 


POLYMER  CONCRETE  MIX  DATA 


Cement  Content 

= 

658  lb/C.''. 

Coarse  Aggregate 

= 

1737  Ib/C.Y. 

(//57  Gravel) 

Absorption 

— 

1.2% 

Specific  Gravity 

— 

2.64 

Fine  Aggregate 

— 

1107  lb/C.Y. 

(natural  sand) 

Absorption 

= 

2.7% 

Specific  Gravity 

= 

2.50 

Water/Cement  Ratio 

= 

.40 

Polymer  Additive 

= 

10  gal./C.Y. 

(manufactured  b 
Dow  Chemical  Co 

Percent  Entrained 

Air 

= 

4.0 

Slump 

= 

3" 

Compressive  Strength 
(28  days) 

= 

4800  psi 

TABLE  2A 


AGGREGATE 


GRADATION 


LIMITS 


#57  Lime  Stone 


Sieve  Size _ %  Passing 


Natural  Sand 


Sieve  Size _ %  Passing 


(2)  2  x  3  x  12  In  (51  x  76  x  305  mm) 

5.3.2  Asphalt  Concrete 

Two  different  sizes  of  crushed  limestone  aggregates  from  the 
local  American  Aggregate  Company  were  utilized  for  the 
manufacture  of  asphalt  concrete  specimens.  The  aggregates  used 
were  #8  crushed  1 imestone  as  coarse  aggregate,  and  1 imestone  dust 
as  fine  aggregate.  The  #8  crushed  aggregate  generally  showed 
sharp,  angular,  and  gritty  particles,  and,  for  the  most  part, 
contained  at  least  one  fractured  face  in  the  particles,  and  were 
reasonably  free  from  excessive  dust  or  other  deleterious 
coatings,  weathered  pieces,  or  excessive  flaky  and/or'  elongated 
pieces.  Measured  water  absorption  of  particle  size  ranged  from 
between  2.9%  and  3.2%.  Aggregate  gradation  conformed  to  FAA 
specif i cat  ions  for  asphalt  concrete  surface  course.  Gradation 
ranges  of  the  #8  crushed  aggregate,  based  on  frequent  samplings, 
are  shown  in  Table  3.  Asphalt  cement  (AC-20)  used  in  this 
investigation  was  obtained  from  Chevron  Asphalt  Company.  Table  2A 
gives  details  about  the  aggregate  gradation  of  #57  limestone  ag¬ 
gregate  and  natural  sand,  used  in  development  of  P-401  mix.  Design 
of  aggregate  blends  for  this  investigations  is  based  on: 


i)  Raw  material  aggregates  and  their  individual  gradings 
11)  General  conformance  with  FAA  specifications  for  P-401  mix 


For  mix  design  and  investigation,  five  levels  of  binder  were 
used,  namely:  5.5,  6,  6.5,  7,  and  7.5%  by  weight  of  total  mix. 

The  optimum  asphalt  content  for  the  selected  aggregate 
gradation  was  determined  using  Marshall  design  procedures  and  75 
blows/face  compaction  efforts.  Optimum  asphalt  content  was 
determined  as  the  average  of  asphalt  content  for  optimum 
stability,  density,  and  4%  air  voids.  Figure  30  shows  the 
Marshal  1  mix  design  properties  for  this  mix. 

5 . 4  Samp l e  Preparat i on 

The  test  specimens  for  testing  were  of  the  following  types: 

(1)  1.5  in  (38  mm)  thick  base  cement  concrete  beam,  notched 

with  teflon  strip  in  the  notch  and  overlaid  with  2  in  (51 
mm)  thick  asphalt  concrete.  The  concrete  beam  had  tack 
coat  with  SS-lh  at  the  rate  of  0.15  gallons  per  squre 
yard.  (0.68  liters  per  square  meter) 

(2)  Same  as  above  with  three  types  of  engineering  fabric  membrane 

placed  on  the  tack  coat,  before  the  ,AC  overlay 


(3)  Same  as  (1)  but  using  2  in  (51  mm)  thick  cement  concrete 


TABLE  3 

GRADATION  OF  P-401  ASPHALTIC  CONCRETE 


2400 


14 


5  5.5  6  6.5  7 

%  AC  by  Wt.  of  Mix 


5  5.5  6  6.5  7 

%  AC  by  Wt.  of  Mix 
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Figure  30.  Marshall  mix  design  data  for  P-401 


(4)  Same  as  (2)  above  but  using  3  In  (76  mm)  thick  overlay 

(5)  Same  as  (1)  and  (3)  above  but  with  3  In  (76  mm)  thick 

over 1 ay 

(6)  Same  as  (1)  but  with  1  in  (75  mm)  A.C.  base  instead  of 

cement  concrete 

(7)  Same  as  (6)  above,  but  with  a  3  in  (76  mm)  overlay 


For  the  fabric  reinforced  beams,  the  lower  part  of  the 
specimen  was  first  manufactured  with  concrete  as  a  beam  of  size 
1. 5x3x2  in  (38x76x610  mm)  and  then  sawed  into  two  parts. 
Tack  coat  with  SS-lh  was  applied  over  the  concrete  beam  and  then 
the  fabric  membrane  was  placed  over  the  sticky  tack  coat 
material.  was  taken  to  maintain  a  1/8  in  (3  mm)  separation 
between  the  saw-rut  pieces  during  fabric  placement. 

After  the  placement  of  fabric  on  the  top  of  tack  coated  base 
beam  the  fabric  surface  was  smoothly  brushed  to  bring  it  into 
complete  contact  with  the  binder.  The  base  beam  with  tack  coat 
and  the  fabric  thus  prepared  was  placed  in  a  beam  mold  and  a 
weighed  quantity  of  hot  mix  poured  over  the  beam.  A  wire  comb 
was  passed  through  the  loose  material  back  and  forth  for  even 
distribution  of  the  mixture  in  the  mold.  The  mixture  in  the  mold 
was  pressed  under  steadily  increasing  load  until  the  asphaltic 
mixture  was  compacted  to  desired  thickness  rather  than  specific 
load  to  insure  that  desired  density  would  be  achieved.  The  mold 
was  then  dismantled  and  the  specimen  was  placed  on  a  stiff 
support,  such  as  a  piece  of  wood  or  steel,  to  await  testing.  All 
precautions  were  taken  to  prevent  bending  and  any  possible  damage 
to  the  beam  sample  prior  to  the  testing.  The  compacted  test 
specimen  was  allowed  to  rool  at  room  temperature  for  a  minimum  of 
24  hours  prior  to  testing. 

5.5  Laboratory  Testing 

5.5.1  Fatigue  Testing 

The  fatigue  experiments  were  conducted  using  a  beam  on 
elastic  foundation  with  geometry  as  shown  in  Figure  28.  The 
selection  of  this  experimental  set-up  was  based  on  a 
two-dimensional  modeling  of  pavement  structure  in  which  a  beam 
representing  the  pavement  is  supported  on  an  elastic  foundation 
representing  the  subgrade.  The  dimensions  of  the  beam  and 
foundation,  as  well  as  the  stiffness  of  foundation,  are  selected 
with  consideration  to  simulating  the  stress  and  strain  at  the 
bottom  of  a  pavement  structure  subjected  to  traffic  loading.  The 
test  set-up  Is  the  same  as  previously  used  by  Resource 
International  Inc.  and  researchers  at  Ohio  State  University  to 
study  the  fatigue  properties  of  asphaltic  mixtures.  The  fatigue 
tests  were  performed  using  a  dynamic  load  function  of  havers ine 
shape.  An  MTS  electro-hydraul ic  testing  system  was  used  to 


generate  the  load  factor.  To  insure  complete  recovery ^0^^ 

sample  before  the  next  load  cycle,  a  res  per  o  '  £  ^  ad 

was  allowed  between  each  load  application.  The  durat Jon  ^  d 
application  in  all  tests  was  kept  **  l  £  and 

Fatieue  test  data  of  samples  tested  at  40  and  72  degre  l 

22  dlg?eef  C)  and  at  different  stress  levels  are  presented  in 

Appendix  A. 

5.5.2  Horizontal  and  Vertical  Stresses 

The  test  set  up  is  shown  schematically  in  Figure  28  and 
pictorially  in  Figure  31.  A  constant  horizontal  pull  was  applie 
LL  Jack  development  In  the  esprit  overlayw.smeaeured 
with  e  micrometer  guage.  Tests  were  conducted  unttlthecrack 

trJSfil  degrees^F  -t  SKS*-*- 

in  Appendix  B. 


CHAPTER  VI 
DATA  ANALYSIS 


6.1  Fatigue  Tests 


The  fatigue  testing  of  the  laboratory  beams  was  conducted 
under  controlled  applied  load  conditions.  Fatigue  analysis, 
however,  requires  the  determination  of  critical  tensile  strain  In 
the  asphaltic  concrete  and  relating  this  strain  to  the  allowable 
number  of  load  applications,  i.e.,  the  fatigue  response  is  given 
by 


-B 

Nf-Ae 


(99) 


where 

Nf  is  the  number  of  load  applications  to  failure 
c  is  the  critical  tensile  strain 

A,B  are  material  constants  that  depend  also  on  temperature 


The  fundamental  assumption  of  mechanistic  stress/strain 
analysis  is  that  Equation  99  can  be  used  to  describe  the  fatigue 
behavior  of  a  particular  material  when  the  critical  tensile 
strain  is  known,  independent  of  how  that  strain  is  developed. 
All  that  is  required  is  to  determine  the  material  constants  (A,  B 
in  Equation  99)  from  some  simple  laboratory  test,  and,  of  course, 
the  critical  strain  in  the  structure  to  be  anallzed. 

The  analysis  model  described  in  Chapter  IV  is  proposed  as 
the  analysis  method  to  determine  the  critical  strains  resulting 
from  loading. 

The  laboratory  test  data  serves  a  two  fold  purpose  in  this 
study: 

(a)  to  determine  the  material  properties  A,  B  in  Equation  99 

(b)  to  verify  the  stress/strain  analysis  model 

The  latter  objective  is  the  reason  for  the  variety  of  conditions 
used  in  fatigue  testing. 

As  was  described  in  the  previous  chapter,  control  beams 


(standard  asphaltic  concrete  overlays)  were  tested  under  the 
fo 1  lowing  cond i t i ons : 


( 1 )  var i ab 1 e  1 oad 

(2)  variable  thickness 

(3)  variable  existing  pavement 

The  significance  of  the  latter  is  that  stress/strain 
concentrations  above  the  crack  are  a  function  of  the  existing 
layer  properties,  i.e.,  the  concentration  is  more  severe  over  a 
cracked  concrete  pavement  than  over  a  cracked  asphaltic  pavement. 

The  test  data  for  the  control  beams  at  72  degrees  F  (22 
degrees  C)  is  shown  in  Figure  32.  The  critical  strains  have  been 
determined  usipg  the  EFRON  program  described  in  Chapter  IV.  As 
can  be  seen  from  this  figure,  the  regression  equation  fits  the 
data  points  very  well.  There  is  some  scatter  in  this  data 
however,  the  correlation  coefficient  (R-square)  for  this  data  is 
0.95  with  a  standard  error  of  estimate  (in  terms  of  log  Nf)  of 
0.125  -  both  values  are  better  than  generally  reported  for 
fatigue  relationships  where  the  only  variable  is  applied  load. 
The  pertinent  fatigue  parameters  are  presented  in  Table  4. 

The  fatigue  data  for  control  beams  at  40  degrees  F  (4 
degrees  C)  is  shown  also  in  Figure  32.  The  variables  for  these 
sets  were  load  and  overlay  thickness;  the  old  pavement  type  was 
kept  constant  as  concrete.  It  can  again  be  seen  that  the 
analysis  model  explains  the  effect  of  the  variables  very 
sat i sfactor i 1 y . 

The  fatigue  data  for  beams  reinforced  with  the  different 
fabrics  are  shown  graphically  in  Figures  33  through  35  and 
summarized  in  Table  4.  The  test  variables  for  these  spec! mans 
were  the  same  as  for  the  control  beams  at  40  degrees  F  (4  degrees 
C),  i.e.,  load  and  thickness  were  varied,  but  all  overlays  were 
over  existing  concrete  pavements. 

It  should  be  noted  that  the  slopes  (the  B  value  of  Equation 
99)  of  the  fabric  reinforced  beams  (indicated  by  dashed  lines  in 
Figures  33  through  35)  are  the  same  as  the  slopes  (B  values) 
obtained  from  regression  analysis  of  the  control  beams.  The 
fatigue  data  for  the  fabric-reinforced  beams  was  analyzed 
assuming  linear  elastic  response,  i.e.,  that  the  stress  and 
strain  are  proportional  to  applied  load  and  that  the  material 
properties  are  independent  of  stress,  as  discussed  in  Chapter  IV. 
Consequently,  the  fabric  must  behave  in  the  same  way  in  the 
system  independent  of  the  applied  load/stress  level. 

Regression  analysis  of  the  fabric  reinforced  fatigue  data 
indicates  that  the  slopes  at  72  degrees  F  (22  degrees  C)  are  not 
drastically  different  from  the  control  beam  slope,  but  are 
substantially  different  at  40  degrees  F  (4  degrees  C). 


TABLE  4 


Laboratory  fatigue  data  for  treatment  B  at  40°F  and  72°F 


Fatigue  Life  (Nf) 

Laboratory  fatigue  data  for  treatment  C  at  4 0°F  and  72°F 


1x10  5xlOJ  lxlO4  5x10^  1x10- 

Fatigue  Life  (Nf ) 

Laboratory  fatigue  data  for  treatment  D  at  40°F  and  ' 


Thus,  the  assumption  of  linear  elasticity  are  not  entirely 
valid;  however,  non-linear  elastic  and/or  elasto-plastic 
analysis  is  beyond  the  scope  of  this  project.  Furthermore,  as  is 
indicated  by  the  standard  error  values  in  Table  4,  the 
assumptions  of  linear  elasticity  do  not  increase  the  standard 
error  values  by  much  -  the  slight  loss  in  precision  is  amply 
compensated  by  the  great  simplification  in  the  analysis  method. 

Of  course,  the  regression  analysis  results  could  be  used  in 
a  phenomenological  model  to  predict  the  fatigue  response.  The 
disadvantage  of  this  approach  is  that  fatigue  response  cannot  be 
predicted  from  fabric  properties  alone  but  would  require 
laboratory  testing.  However,  the  amount  of  testing  is  not 
extensive,  and  the  improved  confidence  level  in  the  predicted 
response  may  outweigh  tne  disadvanatage . 

The  40  and  72  degrees  F  (4  and  22  degrees  C)  test  data  are 
summarized  in  Figures  36  and  37,  respectively.  It  is  appparent 
from  Figure  36  that  fabric  reinforcement  has  a  significant 
beneficial  effect  at  40  degrees  F  (4  degrees  C)  but  that  the 
difference  between  fabrics  is  rather  small.  Figure  37  shows  that 
the  effect  of  fabrics  on  performance  at  72  degrees  F  (22  degrees 
C)  is  not  very  great,  nor  is  the  difference  between  fabrics.  The 
most  probable  explanation  for  this  phenomenon  is  that  although 
the  off-the-shelf  fabric  properties  are  different  among  the 
various  fabrics,  the  fabric-tack  coat  system  properties  are  not 
that  different.  Furthermore,  the  in-situ  fabric-tack  coat  system 
properties  at  72  degrees  F  (22  degrees  C)  are  not  that  different 
from  asphaltic  concrete  properties,  so  that  fabric  effectiveness 
in  relieving  stresses  is  not  great.  At  40  degrees  F  (4  degrees 
C) ,  however,  the  difference  between  fabric  system  and  asphaltic 
concrete  increases,  leading  to  greater  stress  relief  and  fabric 
effectiveness. 

6.2  Simulated  Thermal  Stress  Tests 

The  test  data  using  the  procedure  discussed  in  Section  5.2 
is  presented  in  Appendix  B  in  graphical  form.  The  fundamental 
objective  of  this  test  series  was  to  establish  the  allowable 
joint  opening  at  which  failure  (reflection  cracking)  begins, 
i.e.,  what  is  the  critical  allowable  strain  in  the  asphalt 
overlay,  and  how  various  factors,  such  as  strain  rate,  mean 
temperaure,  overlay  thickness,  reinforcement  type,  affect  this 
value.  The  analysis  of  this  data  is  on  somewhat  more  uncertain 

§  rounds  than  the  fatigue  relationships  since  well-defined  models 
or  prediciting  allowable  critical  strains  are  not  available. 
The  problems  are  complicated  further  by  the  fact  that  since 
asphaltic  concrete  is  a  viscoelastic  material,  loading  rate  is 
likely  to  influence  the  results.  However,  simulation  of  actual 
seasonal  temperature  effects  is  hardly  practical  since  tests 
would  have  to  be  conducted  over  a  several  month  period. 
Therefore,  laboratory  testing  was  conducted  at  loading  rates  that 
are  representative  of  daily  and  weekly  temperature  cycles. 


82 


D  at  72°F 


Fatigue  Life  (Nf) 

Laboratory  fatigue  data  for  control  a,  and  Fabrics  B,  C  &  :D  at  fixed  slope 


Figure  38  is  more  or  less  typical  of  the  creep  data  obtained 
in  this  series  of  tests  and  shows  the  increase  in  joint  opening 
as  a  function  of  time.  This  joint  opening  curve  can  be 
interpreted  to  consist  of  a  linear  (with  time)  and  a  non-linear 
portion,  as  shown  in  the  above  figure.  The  point  at  which  the 
curve  becomes  non-linear  can  be  interpreted  as  the  maximum 
allowable  joint  opening  for  no  cracking.  Small  cracks  were 
observed  on  occasion  at  lower  joint  openings  than  the  allowable 
value  but  for  other  tests  cracks  were  not  observed  until  the 
opening  had  exceeded  the  linear  region,  and  in  general  the  first 
observed  crack  occurred  at  openings  near  the  critical  value 
defined  above. 

Table  5  shows  a  summary  of  the  test  data  for  different 
treatments.  The  data  presented  in  this  table  shows  that  for 
treatment  A  (control),  there  is  very  little  difference  in  the 
allowable  critical  joint  opening  with  thickness  or  with 
temperature.  Figure  39  shows  the  results  of  a  finite  element 
analysis  (using  the  EFRON  program)  of  unreinforced  overlays 
subjected  to  thermally  induced  joint  openings.  As  can  be  seen 
from  this  figure,  the  overlay  modulus  has  an  insignificant  effect 
on  stresses  below  about  200,000  psi  (1.4  MPa),  and  increasing  the 
overlay  thickness  from  2  to  3  in  (51  to  76mm)  decreases  stress  by 
around  7  percent  only.  The  reason  that  critical  stress  is 
insensitive  to  overlay  thickness  is  that  this  stress  is  the 
result  of  stress  concentration  over  the  crack  in  the  concrete 
rather  than  average  stress  in  the  overlay.  Thus  the  measured 
data  for  the  control  beams  agree  quite  well  with  the  theoretical 
analysis . 

The  data  in  Table  5  shows  also  that  reinforcement  of  the 
overlay  with  a  fabric  (Fabrics  B,  C,  D)  substantially  increases 
the  critical  joint  opening,  and  that  the  difference  in 
performance  between  fabrics  is  not  very  significant.  The  same 
conclusion  was  reached  also  from  the  analysis  of  the  fatigue 
data.  The  data  for  fabric  C  at  72  degrees  F  (22  degrees  C)  seems 
to  indicate  superior  performance  of  this  material.  This  data  is 
suspect,  however,  since  the  bond  between  the  concrete  and  the 
overlay  failed  on  many  specimens  before  any  apparent  damage  to 
the  overlay;  the  large  joint  opening  most  likely  represents  some 
slippage  between  the  concrete  and  overlay. 
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Temperature : 


Figure  38.  Typical  joint  opening  creep  curve 


Figure  39.  Maximum  Stress  in  Overlay  Due  to  Joint  Opening 


TABLE  5 

CRITICAL  JOINT  OPENING  FOR  TREATMENT 


Treatment  Temperature  Thickness  Joint  Opening 

F  (  C)  in.  (mm)  mils  (mm) 


A 

40 

(4) 

2 

(51) 

10 

(.25) 

A 

40 

(4) 

3 

(76) 

14 

(.36) 

B 

40 

(4) 

2 

(51) 

23 

(.58) 

B 

40 

(4) 

3 

(76) 

39 

(-99) 

C 

40 

(4) 

2 

(51) 

19 

(.48) 

C 

40 

(4) 

3 

30 

(.76) 

D 

40 

(4) 

2 

20 

(.51) 

D 

40 

(4) 

3 

(76) 

36 

(.91) 

A 

72 

(22) 

2 

(51) 

12 

(.30) 

A 

72 

(22) 

3 

(76) 

11 

(.28) 

B 

72 

(22) 

2 

(51) 

26 

(.66) 

B 

72 

(22) 

3 

(76) 

52 

(1.32) 

C 

72 

l® 

2 

(51) 

85 

(2.16) [ 1  ] 

c 

72 

3 

(76) 

-- 

- [2] 

D 

72 

(22) 

2 

(51) 

30 

(.76) 

D 

72 

(22) 

3 

(76) 

45 

(1.14) 

[1]  Value  erroneous  due  to  slippage  between  concrete  and 
overlay 

[2]  Bond  strength  between  concrete  and  overlay  insufficient 
to  conduct  test 


CHAPTER  VII 


CONCLUSIONS  AND  RECOMMENDATIONS 


7.1  Conclusions 

In  this  research  investigation,  a  mechanistic  methodology  for 
design  of  fabric  reinforced  flexible  overlays  has  been  presented. 
It  has  been  shown  that  the  analytical  model  can  satisfactorily 
predict  the  performance  of  flexible  overlays  of  rigid  airport 
pavements.  To  validate  this  model,  laboratory  fatigue  tests  have 
been  carried  out  under  controlled  load  and  temperature 
environments.  The  three  engineering  fabrics  used  in  this 
laboratory  investigation  correspond  to  low,  medium  and  high 
tensile  moduli  fabrics.  The  following  major  conclusions  are  drawn 
from  this  study. 

1.  The  results  of  the  laboratory  fatigue  tests  clearly 
demonstrates  the  effectiveness  of  fabrics  in  retarding  the 
formation  of  reflective  cracking  of  flexible  overlays  of 
rigid  airport  pavements . 

2.  The  beneficial  effects  of  fabrics  are  greatest  at  the 
lower  temperatures  where  most  help  is  needed.  At  these 
temperature  levels,  it  has  been  shown  that  the  allowable 
strain  values  are  lower,  requiring  a  greater  contribution 
by  fabric.  The  results  of  the  temperature  simulation 
tests  also  demonstrate  the  beneficial  effect  of  using 
fabric  reinforcement  in  retarding  the  formation  of 
reflective  cracks. 

3.  For  the  range  of  variables  investigated,  the  test  results 
show  very  little  difference  in  the  performance  of  the 
three  engineering  fabrics  used  in  this  study.  However, 
the  results  of  this  research  study  is  not  applicable  to 
the  metallic  reinforcement  and  other  woven  fabric  types. 

4.  In  this  study,  an  SS-1H  emulsion  was  used  as  a  tack  coat 
with  application  temperature  of  77  F  and  application  rate 
of  0.15  gallon  per  square  yard  (0.68  liter  per  square 
meter).  The  selection  of  optimum  quantity  of  tack  coat 
and  tack  coat  type  shall  depend  on  the  Concrete  Surface 
Conditions  as  well  as  the  asphalt  retention  of  the 
fabrics . 

5.  An  optimum  fabric  performance  can  only  be  achieved  under 
limited  horizontal  and  vertical  displacement  of  underlying 
concrete  pavement.  The  horizontal  and  vertical 
displacements  under  load  and  temperature  are  known  to 
cause  reflective  cracking  of  flexible  overlays.  It  is 
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shown  that  at  the  40  F,  the  measured  values  of  the  joint 
opening  above  which  a  crack  will  develop  ranges  between  10 
mils  to  14  mils  (0.25  to  036  mm)  for  a  no  fabric  overlay, 
and  ranges  between  19  mils  to  36  mils  (0.48  to  0.91  mm) 
for  the  fabric  reinforced  overlays.  Further,  at  72  F 
(22C)  a  threshhold  value  of  joint  opening  is  about  11  mils 
(0.28  mm)  for  a  no  fabric  overlay  while  it  ranges  from  26 
mils  to  52  mils  (0.66  to  1.32  mm)  for  the  fabric 
reinforced  overlays. 

6.  As  a  general  conclusion,  the  reinforcement  of  the  overlay 

with  a  fabric  substantially  increases  the  critical  joint 
opening,  and  that  the  difference  in  performance  between 
fabrics  tested  in  this  study,  is  not  very  significant.  On 
the  other  hand  the  large  joint  opening  most  likely 
represents  some  slippage  between  the  concrete  and  overlay 
based  on  the  laboratory  observation  that  the  bond  between 
the  concrete  and  the  overlay  failed  on  many  of  these 
specimens  before  any  apparent  damage  occurred  to  the 

overlay. 

7.  For  reinforced  overlays  subject  to  thermally  induced  joint 

openings,  a  threshhold  modulus  value  of  200,000  psi  (1.4 
mpa)  was  found  using  EFRON  program,  below  which  the 

modulus  has  no  effect  on  the  maximum  stress  in  the 

overlays . 

In  summary,  it  should  be  pointed  out  that  the  above  conclusions 
could  be  slightly  extrapolated  beyond  laboratory  conditions  used 
in  this  study.  The  comparison  of  the  test  results  with  the 

analysis  using  Mechanistic  Model  EFRON  program  shows  that  the 

analysis  techniques  utilized  in  that  program  realistically  model 
the  behavior  of  the  overlay  system  in  the  laboratory  tests,  and 
there  is  every  expectation  that  the  model  (EFRON)  will 

satisfactorily  predict  the  performance  of  overlays  of  actual 
runways  also. 

The  verification  of  the  applicability  of  the  EFRON  program  to 

runway  overlays  is  the  major  task  of  Volume  II  of  this  study. 

EFRON  program  is  available  at  the  FAA  library. 


7 . 2  Recommendations 


Based  on  this  research  study,  it  is  recommended  that; 

1.  The  type  of  fabric  membrane  to  be  used  to  improve  the 
fatigue  life  of  pavement  shall  be  non-woven  having  a 
tensile  strength  of  not  less  than  90  lbs.  (41  kgs) 
determined  according  to  ASTM  D  1682  and  a  density  in 
the  range  of  3  to  5.5  ozs.  per  square  yard  (40  to  70 
gms  per  square  yard). 

2.  Non-woven  fabric  membranes  are  to  be  used  where  the 
horizontal  displacement  is  not  in  excess  of  50  mils 
and  vertical  displacement  is  not  in  excess  of  20  mils. 

3.  Tack  coat  to  be  used  for  most  non-woven  fabrics  shall  be 
emulsified  asphalt  with  the  rate  of  application  ranging 
from  0.15  to  0.30  gallons  per  square  yard.  The  optimum 
tack  coat  quantity  shall  depend  on  the  type  of  fabric 
membrane  used  and  the  surface  condition  of  the  airport 
pavement  to  be  overlaid. 

4.  The  asphalt  concrete  overlay  thickness  when  fabric 
membranes  are  used  shall  not  be  less  than  3  inches 
(75  mm). 

Construction  specifications  to  be  adopted  for  fabric  reinforced 
overlays  will  be  detailed  in  Volume  II  of  this  report  and  the 
above  recommendations  may  be  subject  to  change  based  on  the 
final  findings  of  this  research  study. 


APPENDIX  A 


LABORATORY  FATIGUE  TEST  DATA  ON  BEAMS 

All  beams  3  In  (76  mm)  wide  x  24  In  (610  mm)  long 

Treatment  A  -  control,  or  conventional  overlay 

Treatment  B  -  overlay  reinforced  with  fabric  B 

Treatment  C  -  overlay  reinforced  with  fabric  C 

Treatment  D  -  overlay  reinforced  with  fabric  D 

Note:  For  properties  of  fabrics 

refer  to  Table  1:  page  67 


TABLE  6 


FATIGUE  TEST  DATA  AT  72  DEGREES  F 


V 

Beam  Base 

Base 

Overlay 

Specific 

Number  Thickness 

Type 

Thickness 

Gravity  ] 

I 

i 

(inches) 

(inches) 

1 

2 

3 

4 

5 

A-2-6 

7 

8 
9 

A-2-10 
A-2-11 
A-2-12 
1 

2 
3 

5 
5 
7 
3 
3 

B-2-10 
B-2-11 
B-2-I2 
1 

2 


D-2-2 


cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
cone 
C 
C 


TABLE  6  (continued) 

FATIGUE  TEST  DATA  AT  72  DEGREES  F 


Beam 

Number 

Base 

Thickness 

(inches) 

Base 

Type 

Overlay 

Thickness 

(inches) 

Specific 

Gravity 

Test 

Load 

(lbs) 

Fatigue 

Life 

CB-15 

1.01 

A.C . 

2.00 

2.36 

150 

46.00 

CB-1 

1.0 

A.C . 

2.00 

2.36 

200 

4.40 

CB-2 

1.05 

A.C. 

2.00 

2.36 

200 

11.40 

CB-3 

1.0 

A.C. 

2.00 

2.35 

200 

8.79 

CB-4 

1.0 

A.C. 

2.00 

2.35 

200 

2.70 

CB-5 

1.0 

A.C. 

2.00 

2.36 

200 

1 .85 

CB-6 

1.0 

A.C. 

2.00 

2.36 

200 

1.60 

CB-10 

1.0 

A.C. 

2.00 

2.37 

250 

1 .02 

CB-1 1 

1.1 

A.C. 

2.00 

2.37 

250 

1 .09 

CB-12 

1.1 

A.C. 

2.05 

2.36 

250 

0.96 

CB-16 

1.1 

A.C. 

2.05 

2.35 

250 

2.50 

CB-17 

1.0 

A.C. 

2.00 

2.33 

250 

1.60 

CB-1 8 

1.05 

A.C. 

2.00 

2.32 

250 

1.20 

A-3-1 

1.5 

cone . 

3.00 

2.37 

150 

48.00 

A-3-2 

1.5 

cone . 

2.98 

2.37 

150 

44.00 

A- 3- 3 

1.5 

cone. 

2.97 

2.37 

200 

12.80 

A-3-4 

1.5 

cone. 

2.99 

2.37 

200 

14.30 

A-3-5 

1.5 

cone . 

3.00 

2.38 

250 

3.70 

A-3-6 

1.5 

cone . 

3.01 

2.35 

250 

3.60 

B-3-1 

1.5 

cone. 

3.00 

2.36 

150 

137.8 

B-3-2 

1.5 

cone . 

3.00 

2.36 

150 

251 .4 

B-3-3 

1.5 

cone . 

3.00 

— 

200 

122.0 

B-3-4 

1.5 

cone . 

3.00 

— 

200 

112.0 

B-3-5 

1.5 

cone . 

3.00 

2.37 

250 

26.5 

B-3-6 

1.5 

cone . 

3.00 

2.35 

250 

18.8 

C-3-1 

1.5 

cone . 

3.00 

2.36 

150 

42.5 

C-3-2 

1.5 

cone . 

3.00 

2.38 

150 

58.0 

C-3-3 

1.5 

cone . 

2.99 

2.38 

200 

44.0 

C-3-4 

1.5 

cone . 

3.01 

2.37 

200 

60.4 

C-3-5 

1.5 

cone . 

3.03 

2.36 

250 

18.2 

C-3-6 

1.5 

cone . 

3.01 

2  .35 

250 

14.0 

D-3-1 

1.5 

cone . 

3.00 

2.36 

150 

301  .5 

D-3-2 

1.5 

cone . 

3.00 

2.37 

150 

326.0 

D-3-3 

1.5 

cone . 

3.00 

2.37 

200 

110.3 

D-3-4 

1.5 

cone . 

3.00 

2.32 

200 

140.5 

D-3-5 

1.5 

cone . 

3.00 

2.36 

250 

27.9 

D-3-6 

1.5 

cone . 

3.00 

2.37 

250 

28.2 

TABLE  7 


FATIGUE  TEST  DATA  AT  40  DEGREES  F 


Beam  Base  Base  Overlay  Specific 

Number  Thickness  Type  Thickness  Gravity 

(inches)  (inches) 


Test  Fatigue 
Load  Life 
(lbs) 


A-2-13 

2.0 

A-2-14 

2.0 

A-2-15 

2.0 

A-2-16 

2.0 

A-2-17 

2.0 

A-2-18 

2.0 

A-2-19 

2.0 

A-2-20 

2.0 

A-2-21 

2.0 

B-2-13 

2.0 

B-2-14 

2.0 

B-2-15 

2.0 

B-2-16 

2.0 

B-2-17 

2.0 

B-2-18 

2.0 

B-2-19 

2.0 

B-2-20 

2.0 

B-2-21 

2.0 

C-2-7 

2.0 

C-2-8 

2.0 

C-2-10 

2.0 

C-2-11 

2.0 

C-2-13 

2.0 

C-2-14 

2.0 

C-2-15 

2.0 

D-2-10 

2.0 

D-2-11 

2.0 

D-2-12 

2.0 

D-2-7 

2.0 

D-2-8 

2.0 

D-2-9 

2.0 

D-2-13 

2.0 

D-2-14 

2.0 

A- 3- 7 

1.5 

A-3-8 

1.5 

A-3-10 

1.5 

A-3-11 

1.5 

A-3-13 

1.5 

A-3-14 

1.5 

B-3-7 

1.5 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 

cone 


2.0 

2.36 

150 

10.2 

2.0 

2.32 

150 

11.1 

2.05 

2.32 

150 

10.9 

2.0 

2.37 

200 

6.0 

108.34 


365.90 


Test  Temperature:  40Op  (4°4  C) 
Constant  Load:  140  lbs  (621. 6N) 


ature:  40°F  {4.4 

ad:  155  lbs.  (688 
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Test  Temperature:  40°F  (4.4°C) 
Constant  Load:  140  lbs.  {621. 6N) 
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Test  Temperature:  40°F  (4.4°C) 
Constant  Load:  155  lbs.  (688. 2N) 
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ELAPSED  TIME  (HRS) 


Test  Temperature:  72°F  (22°C) 
Constant  Load:  5  2  lbs.  (230. 9 N) 


ELAPSED  TIME  (HRS) 


72°F  (22°C) 


ELAPSED  TIME  (HRS) 
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Note 


Test  Temperatures  40°F  (4.4°C) 
Constant  Load:  155  lbs.  (688. 2N) 
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APPENDIX  C 

MECHANISTIC  METHODOLOGY  FOR  AIRPORT 
PAVEMENT  DESIGN  WITH  ENGINEERING  FABRICS 

PHASE  II:  USER  GUIDE 
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EFRON  USERS  MANUAL 


1  .  PROGRAM  CAPAC I TY 


EFRON  is  a  computer  program  designed  to  analyze 
stress-strain  distribution  for  existing  pavements  covered  with 
rigid  or  flexible  overlays.  It  is  based  on  the  theory  described 
in  Volume  1  of  this  report. 

An  existing  pavement  structure  can  be  supported  by  an 
N-layered  elastic  foundation  (N<3)  which  is  represented  in  this 
program  by  a  stiffness  matrix  computed  based  upon  Burmeister's 
solution  of  3  dimentional  semi -infinite  space.  Overlay  structure 
can  contain  an  engineering  fabric  which  can  act  as  a  bond  breaker 
and/or  as  a  reinforcement.  EFRON  utilizes  a  two-d i ment i ona 1 
finite  element  approach  and  has  the  following  capabilities: 

1.  Slabs  can  have  a  variable  length  and  can  be 
divided  for  regular  or  irregular  infinite 
element  mesh,  which  is  automatically  created 
by  the  program. 

2.  Boundary  conditions  are  assumed  by  the  program 
depending  on  assumed  load  type,  support 
structure  and  boundary  condition  option. 

3.  Three  types  of  the  foundations  can  be 
specified  as  an  option: 

a)  solid  elastic 

b)  Winkler 

c )  Rigid 

4.  Three  load  types,  or  any  combination  of  these, 
are  avai 1 abl e : 

a)  Vertical  Load  (mechanical) 

b)  Vertical  temperature  gradient 
(daily  temperature  change) 

c)  Uniform  temperature  load 
(seasonal  temperature  change) 

5.  Number  of  support  layers  can  be  varied  from 
1  to  3 . 


6. 


Iteration  procedure  can  be  specified  as  an 
option  due  to  partial  contact  and/or  elasto- 
plastic  shear  resistance  between  slab  and 
foundation.  A  special  bond  type  element  has 
been  applied  for  this  feature. 

7.  Engineering  fabric  is  simulated  by  N- layered 
stripe  which  can  be  placed  at  any  position  in 
the  overlay.  The  number  of  fabric  layers  (N) 
can  vary  from  0  (no  fabric)  to  4.  Each  layer 
of  this  stripe  can  have  separate  material  and 
geometric  properties. 

8.  All  material  properties  of  the  structure  may  be 
orthotropic  and/or  temperature  dependent.  This 
feature  is  extremely  important  for  fabric  whose 
modulus  in  vertical  directon  is  different  from 
hor i zonta 1 . 

9.  Transverse  joint  may  be  plain  or  doweled. 

10.  Dowel  bar  looseness  can  be  specified  between 
0  and  8  mi  1  (0  and  .2  mm)  by x a  stiffness 
reduction  function  developed  by  Resource 
International  Inc. 

11.  Variable  size  void  under  the  slabs  can  be 
specified. 

12.  Plane  stress  or  plane  strain  analysis  can  be 
spec  if led. 

13.  New  pavement  or  overlay  design  can  be  performed. 

14.  Up  to  10  different  cracks  can  be  specified  in 
different  locations  of  the  slabs. 

A  Simplified  flow  chart  of  the  EFRON  computer  program  is 
shown  in  F i gure  1 . 

2.  FINITE  ELEMENTS  MESH  AND  BOUNDARY  CONDITIONS 


A  finite  element  mesh,  as  well  as  boundary  conditions,  will 
be  generated  automatically  by  the  program  depending  on  imput 
parameters  (see  below  for  input  instructions  and  mesh  examples). 
However,  the  following  recommendations  for  the  use  of  boundary 
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Compute  residual  load  vector 
Pi 


Resolve  Kfr  =  p 


- Are  convergence 
criteria  satisfied? 


Figure  1  (cont'd).  Flow  Chart  of  EFRON 


conditions  option  and  mesh  type  options  should  be  considered: 

1.  Use  regular  mesh  in  case  of  no  mechanical  load 
applied  on  the  structure  (IR=0). 

2.  Restrict  horizontal  movement  on  each  edge  in  the 
above  case  (IBC=0).  This  boundary  condition  means 
that  this  is  a  double  symmetry  case. 

3.  Using  a  partial  restriction  in  x,y  directions  is 
recommended  in  case  of  a  simulation  of  the  support 
of  the  infinite  or  very  long  slab.  This  is  a 
Winkler  type  of  spring  and  its  axial  stiffness 
calculation  is  not  very  precise. 

4.  Do  not  use  a  very  small  or  very  large  element  size 
in  case  of  regular  mesh.  A  very  small  size 
significantly  increases  computation  time.  A  very 
large  size  decreases  result  credibility.  Recommended 
size  is  1/12  -  1/15  of  a  single  slab  1 ength . 

5.  Assume  a  minimum  vertical  element  size  (DLS1ZE) 
equal 

DLSIZE  *  ,2*TH( 1 ) 

where 

TH ( 1 )  overlay  thickness 

This  assumption  will  result  in  a  more  correct 
stress  above  the  joint. 

Figure  2  shows  a  two-d 5 ment i ona 1  representation  of  the 
overlayed  pavement  structure  used  in  the  finite  element  (EFRON) 
program.  This  figure  shows  also  the  definition  of  several  of  the 
input  variables.  Figure  3  shows  a  possible  FEM  mesh  for  this 
pavement  which  is  an  irregular  mesh  in  this  case.  As  was  stated 
above,  the  irregular  mesh  is  recommended  in  the  case  of  wheel 
load  (mechanical  load),  but  a  regular  mesh  should  be  used  with 
thermal  loads. 

The  information  required  as  input  is  described  in  Table  1  in 
conjunction  with  Figures  2  through  4.  The  makeup  of  the  data 
deck  is  shown  in  Figure  5.  It  may  seem  from  these  figures  as  of 
a  great  deal  of  data  is  required;  however,  the  imput  data  is  in 


I  FT ION  OF  INPUT  DAT 


16-20  Number  of  Iterations. 

0  -  No  iterations. 


LSHOUL  =  in  Fig.  2. 

1  -  Shoulder .  if  l  should  = 

IBC  must  be  1 


Table  11 
cont'd 

DESCRIPTION  OF  INPUT  DATA 


Table  *1 
cont'd 

DESCRIPTION  OF  INPUT  DATA 


cont'd 

DESCRIPTION  OF  INPUT  DATA 


2  ELSIZE  11-20  Minimum  Element  Size.  (in. 


CARD  FORMAT  VARIABLE  COL.  DESCRIPTION  &  VALUE  RANGE  COMMENTS 


5  E(I, 4)  41-50  Shear  Modulus  of  Fabric  Layer  (I) 


CARD  FORMAT  VARIABLE  COL.  DESCRIPTION  6,  VALUE  RANGE  COMMENTS 


3  V(9,l)  21-30  Subbase  Poisson's  Ratio 


5F10.0  1  SNORM  1-10  Normal  Stiffness  for  Joint 

Element  (psi) 

2  SSHR  11-20  Shear  Stiffness  for  Joint  Elemeni 

(psi) 


Table  *1 
cont'd 

DESCRIPTION  OF  INPUT  DATA 


615  1  NELM  1-5  Beginning  element  number  (see  Use  this  Card 

mesh) .  only  if  10  =  ! 


cont ‘d 

DESCRIPTION  OP  INPUT  DATA 


80  COMMENTS 


general  not  difficult  to  generate  once  the  user  has  decided  on 
the  appropriate  options. 


3.  RECOMMENDATIONS  OF  INPUT  MATERIAL  PROPERTIES 


As  was  stated  previously,  the  slab  support  layers  are 
characterized  by  using  elastic  layer  theory;  therefore,  elastic 
material  properties  are  required  (Young's  modulus,  Poisson's 
ratio,  and  layer  thickness)  for  each  foundation  layer  as  well  as 
for  the  concrete  slab.  Additionally,  flexural  strength  (28  day, 
3  point  loading)  is  required  for  concrete. 


(a)  The  flexural  strength  (modulus  of  rupture) 
of  most  airport  concretes  is  around  600  - 
800  psi  (4.1  to  5.5  MPa).  This  value  is 
generally  available  from  test  data  but 
may  be  estimated  from  unconfined  compressive 
strength  tests  using  the  following  ACI 
relationship: 


where 


fc  =  a  f'c 


( 1) 


fc  is  the  flexural  strength 

f'c  is  the  compressive  strength  of  standard 

6  X  12  in.  (152  X  305  mm)  cylinders 

a  i s  a  constant  rang i ng  from  7  to  1 0  for 

British  units  and  0.58  to  0.83  for 

International  units. 

The  lower  value  should  be  used  with  high- 
strength  concrete  and  the  higher  value  for 
low  strength  concrete. 

(b)  The  dynamic  (tangent)  modulus  of  highway 

concretes  is  around  5  million  psi  (34.5  GPa) 
within  a  range  of  3.5  to  6  million  psi 
(24.1  to  41.4  GPa).  In  the  absence  of  specific 
test  data,  it  may  be  estimated  from  the  following 
ACI  relationship: 


Ec  *  43 


*f 'c  (in  psi ) 


where 


(2) 


1.5 


Ec  *  0.056 


f'c  (in  kPa) 


Ec  is  the  tangent  modulus 

is  the  concrete  unit  weight 
f'c  is  the  compressive  strength 


Equation  2  is  generally  given  for  the  secant  modulus.  Since 
the  tangent  modulus  is  usually  20  to  30  percent  higher  than  the 
secant  modulus,  the  coefficients  in  equation  2  have  been 
increased  by  30  percent  over  the  values  recommended  by  AC  I  for 
secant  modulus  (for  a  more  conservative  design)  because  critical 
stresses  in  concrete  increase  with  increasing  concrete  modulus. 

(c)  The  Poisson's  ratio  values  in  Table  2  are 
recommended .  The  ana lysis  is  rather  i nsens i t i ve 
to  Poisson's  ratio  so  that  precise  values  are 
not  required. 

(d)  Subgrade  modulus  may  be  obtained  from 
laboratory  tests  in  tri axial  compression, 
or  may  be  est i mated  from  A ASHTO  so i 1 
classification,  soil  support  values,  or 
modulus  of  subgrade  reaction,  using 
Figure  6  to  estimate  the  CBR  value,  and 
Equation  3  to  compute  the  modulus. 


Es  = 

1500 

CBR 

( in 

psi ) 

Es  = 

10.3 

CBR 

( in 

MPa) 

(3) 

(e)  Granular  base/subbase  moduli  may  be  obtained 
from  laboratory  tests  in  triaxial  compression, 
or  may  be  estimated  from  Equation  4s 

0.837 

En-1  =  11.06  EN  (In  psi )  (4) 


147 


*. 


/,  s 


TABLE  2 

RECOMMENDED  POISSON'S  RATIOS 

Material  Poisson's  Ratio 


Concrete 

0.15 

Asphaltic  concrete 

0.40 

Granular  base 

0.37 

Cement-treated  base 

0.20 

Lime-stabilized  soil 

0.35 

Subgrade 

0.45 

CALIFORNIA  BEARING  RATIO  (CBR) 


v 


.SOIL  SUPPORT  VALUE  (SSV) 


Figure  6.  Relation  between  soil  properties  and  CB12 


where 


EN  a  modulus  of  the  nth  layer  above  it. 

EN-1  =  modulus  of  the  layer 

This  equation  is  a  compromise  between  the  Shell 
COE,  and  Kentucky  models  and  recognizes  the  fact 
that  the  degree  of  compaction  of  granular  materials 
depends  on  the  modulus  of  the  underlying  layer.  If 
both  a  subbase  and  a  base  layer  is  used.  Equation  4 
should  be  applied  twice,  first  to  determine  the 
subbase  modulus  from  subgrade,  then  the  base  modulus 
from  subbase  modulus.  Sensitivity  analyses  using 
both  elastic  layer  theory  as  well  as  the  RISC  program 
show  that  the  critical  stress  in  concrete  is  only 
slightly  dependent  of  the  value  of  the  base/subbase 
modulus;  therefore,  exact  values  for  these  layers 
are  not  required. 

(f)  The  modulus  of  asphalt  bases  is  temperature 
dependent  and  should  be  estimated  at  mean  annual 
temperature.  This  mean  temperature  depends  on  the 
slab  thickness  and  the  base  thickness  as  well  as 
on  the  mean  annual  air  temperature,  but  using  mean 
annual  air  temperature  will  be  adequate  in  most 
cases.  Figure  7  may  be  used  in  absence  of  test 
data . 

(g)  The  moduli  of  cement-treated  bases  are  quite 
variable,  depending  on  aggregate  type  and  the 
cement  contents;  therefore,  the  use  of  laboratory 
test  data  is  recommended. 


The  above  discussion  should  provide  the  user  some  guidelines 
in  the  selecton  of  imput  parameters  required  for  this  design 
program.  However,  since  the  program  has  the  capability  of 
analyzing  a  wide  variety  of  programs,  some  planning  and 
forethought  is  recommended. 


rtuPCRATUAC  ,  •  F 


Figure  7.  Assumed  temperature  dependence  of  Young’ 
modulus  of  AC  pavements  and  AC  base  materials 
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